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ABSTRACT 

Energy  absorbing  mechanisms  in  sand  subjected  to  one-dimensional 
compression  are  reported.  For  pressures  below  crushing  of  the  grains, 
a  granular  medium  of  equal  radii  elastic  spheres  in  a  face-centered  cubic 
array  is  analyzed.  Expressions  are  obtained  for  the  axial  stress- strain 
curve,  constrained  modulus,  coefficient  of  earth  pressure  at  rest,  and 
relationship  between  absorbed  energy  and  input  energy  for  one  cycle  of 
loading.  The  energy  absorbed  as  a  result  of  crushing  is  considered  by 
analyzing  statistical  relationships  between  changes  in  grain  size  distribution 
curves  and  the  new  surface  areas  created.  An  apparatus  is  described  which 
has  the  capability  of  maintaining  conditions  at  zero  radial  strain  under  in¬ 
creasing  axial  stress.  The  lateral  stresses  developed  under  these  conditions 
are  measured.  Preliminary  experimental  results  are  presented  for  one 
sand  which  show  the  variation  of  the  coefficient  of  earth  pressure  at  rest 
and  the  stress- strain  relationships  with  initial  void  ratio,  overconsolidation 
ratio,  and  strain  rate.  A  correlation  of  theory  and  test  results  is  presented. 
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CHAPTER  1 

urraoDocTiOH 

1.1  Object  and  Scope 

The  object  of  this  study  Is  to  define  the  mechanisms  by  vhlch 
granular  materials  such  as  dry  sands  and  gravels  absorb  energy  vhen  sub¬ 
jected  to  an  applied  state  of  stress.  The  study  therefore  Is  an  Investi¬ 
gation  of  the  stress-strain  behavior  and  hysteresis  effects  In  cohesion¬ 
less  granular  materials  subjected  to  Increments  of  stress  extending  Into 
the  higher  pressure  ranges. 

The  usual  concept  of  a  stress -strain  relationship  for  elastic 
solids  does  not  apply  to  granular  materials.  The  Ideeil  group  of  materials 
classified  as  elastic  solids  exhibit  a  linear  scress-strain  relation.  In 
such  a  material  the  stress-strain  properties  may  be  described  by  elastic 
constants.  A  grantilar  medium  such  as  soil  Is  entirely  different  how¬ 
ever^  in  that  constants  cannot  be  used  to  describe  the  stress -strain 
properties  of  the  medium  even  If  the  material  properties  of  the  p€urtl- 
cles  and  the  packing  arrangement  are  knovn.  The  stiffness  of  a  granular 
medlisn  Is  not  only  dependent  upon  the  material  composition  of  the  peurtl- 
cles  and  the  packing  arrangement  (relative  density )>  but  the  stiffness 
Is  also  a  non-linear  function  of  the  stress  tensor.  The  dependence  of 
the  stress -strain  relation  can  be  Illustrated  by  considering  various 
states  of  stress  on  the  cylindrical  sample  of  soil  shovn  In  Fig.  1.1. 

Curve  1  shovs  an  axial  stress-strain  cxuve  for  a  sample  of  soil  idilch 
vas  compressed  hydrostatically  with  c  »  a  .  The  concave  upward  stress- 
strain  relationship  Is  a  typical  relationship  because  as  the  soil  be¬ 
comes  denser  the  resistance  to  volune  change  Increases.  Curve  3  illus¬ 
trates  an  axial  stress -strain  curve  for  a  sample  under  constant  radial 


stress  ^ich  Is  deformed  ty  increasing  the  axled  stress.  Ihe  resulting 
stress-strain  curve  is  concave  dovmward  indicating  that  the  resistance 
to  defoimatlon  Increases  at  a  decreasing  rate.  This  behavior  is  typical 
vrtiere  the  primary'  resistance  to  deformation  is  a  sheeurlng  resistance 
rather  than  a  resistance  to  volume  change.  Chirve  2  Is  a  one  dimensional 
compression  curve  ^ere  the  axial  stress  is  increased  under  the  condi¬ 
tions  of  zero  radial  strain.  In  this  situation  the  lateral  stresses  are 
not  controlled  and  are  statically  Indeterminate.  Uie  concave  upward 
stress -strain  curve  shows  that  the  resistance  to  deformation  in  one 
dimensional  compression  is  primarily  a  resistance  to  a  change  in  volume. 
These  examples  Illustrate  the  various  types  of  behavior  manifested  by 
soil  and  demonstrates  the  dependence  of  the  stress-strain  relation  on 
the  applied  state  of  stress.  The  general  stress -strain  relations  in  a 
granular  material  are  therefore  very  complicated,  especially  %rtien  the 
shearing  stresses  predominate  as  in  curve  3  and  particles  rotate,  slide, 
and  a  rearrangement  of  peirtlcles  is  constantly  taking  place. 

In  order  to  eliminate  the  variables  Introduced  by  gross  re¬ 
arrangement  and  other  uncertainties  which  occ\ir  under  large  shearing 
stresses,  this  study  was  limited  to  those  cases  idiere  the  shearing 
stresses  are  small  with  respect  to  the  normal  stresses.  Die  behavior 
in  one  dimensional  compression  has  been  emphasized  because  this  case  is 
of  practical  Importance  to  the  Air  Force.  Ihe  one  dimensional  case  is 
a  reasonable  representation  for  an  istportant  group  of  protective  construc¬ 
tion  problems;  those  ^ere  the  air  blast  induced  ground  shock  is  of  pri¬ 
mary  Importance.  In  this  condition  of  loading  the  radial  extent  of  the 
loaded  area  is  large  in  comparison  to  the  thickness  of  highly  compressible 
surface  soil  and  near  surface  unconsolidated  rock  strata.  Under  this 
condition  of  loading,  the  strata  are  latercdly  confined  such  that  the 
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only  displacement  that  can  occur  Is  In  the  direction  of  stress  wave  pro* 
pagatlon. 

A  typical  stress -strain  curve  for  a  granular  material  such  as 
sand  subjected  to  static  loads  In  one  dimensional  compression  Is  shown 
qualitatively  in  Fig.  1.2.  This  curve  shows  that  the  stress-strain  be¬ 
havior  and  consequently  the  energy  absorption  mechanisms  are  dependent 
to  a  large  extent  on  the  magnitude  of  pressure. 

The  behavior  In  the  very  low  stress  ranges,  (Region  1,  Fig. 
1.2),  reflects  a  rearrangement  of  the  grains  and  the  stress-strain  curve 
Is  concave  downward.  Ihe  energy  absorbed  In  the  compaction  process  Is 
of  course  non -re cover able;  however,  if  the  pressure  subsequently  becomes 
quite  large  this  energy  loss  Is  quite  small  compared  to  losses  caused 
by  other  phenomena. 

As  the  load  Is  Increased  the  particles  begin  to  lock  together 
and  become  a  stable  matrix  of  elastic  particles  (Region  2).  In  this 
higher  stress  region  the  strains  that  take  plMe  in  the  material  are  due 
primarily  to  the  deformation  of  the  particles  at  the  points  of  contact. 
There  is  of  course  some  rearrangement  continuing  In  this  region  but 
studies  Indicate  that  the  rearrengement  Is  not  significant  at  this  pres¬ 
sure  level.  The  behavior  in  this  region  Is  essentially  a  non-linear 
elastic  behavior  In  that  there  Is  little  permanent  set  resulting  from 
an  Increment  of  stress  applied  and  removed  within  Region  2.  This  be¬ 
havior  can  be  analytically  studied  by  means  of  the  Hertz-Mlndlln  con¬ 
tact  theory  which  Is  discussed  In  detail  In  Chapter  3*  A  study  of  this 
type  was  made  on  a  face  centered  array  of  equl-radll  spheres.  A  theo¬ 
retical  stress-strain  curve  was  derived  for  the  array  of  si^erea  In 
hydrostatic  compression  and  the  concepts  were  extended  to  obtain  a 


stress -strain  relationship  for  the  array  subjected  to  the  bouzidary  condi¬ 
tions  concomitant  with  one  dimensional  compression.  It  was  found  that  in 
one  dimensional  compression  it  is  possible  to  have  energy  absorbed  due  to 
friction  at  the  contacts  tdien  no  permanent  strain  results  from  a  cycle  of 
loading.  Ihe  question  arises  as  to  vhether  such  a  theory  based  on  the 
theories  of  contact  stresses  is  applicable  to  soils  because  of  the  magni¬ 
tude  of  the  stresses  and  the  small  radii  of  the  grains.  Appendix  A  treats 
this  topic  and  shovs  that  for  the  stress  ranges  of  interest^  namely  those 
in  Region  2  prior  to  particle  crushing^  the  contact  theory  is  applicable. 
For  certain  cases  other  than  complete  confinement  vlth  various  amotints  of 
lateral  deformation  permitted^  it  was  necessary  to  obtain  nimierlcal  solu¬ 
tions  to  the  theory  of  granular  media  since  a  closed  form  solution  was  not 
possible.  Appendix  B  outlines  the  work  carried  out  and  the  program  coded 
for  a  digital  computer. 

As  the  stress  continues  to  increase  to  the  neighborhood  of  one 
to  five  thousand  psl  the  contact  forces  become  so  large  that  the  particles 
begin  to  crush.  The  crushing  is  accompanied  by  a  certain  amount  of  re¬ 
arrangement  which  is  reflected  in  the  concave  downward  stress-strain  curve 
shown  in  Region  3*  Energy  is  absorbed  in  this  rajoge  by  the  creation  of 
new  surface  area  during  crushing  and  energy  is  dissipated  as  beat  due  to 
rearrangement  of  the  new  particles  into  a  denser  configuration.  Both  of 
these  mechanisms  absorb  energy  idilch  is  no  longer  in  the  form  of  recover¬ 
able  strain  energy. 

As  the  pressure  continues  to  increase,  the  particles  lock  again 
and  the  curve  again  becomes  concave  upward  as  shown  in  Region  4,  Fig.  1.2. 
In  general  the  stress -strain  curve  in  Region  4  tends  to  follow  the  same 
general  pattern  ^Ich  took  place  in  Region  2  before  crushing  except  that 
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the  particles  are  smaller  and  more  angular  than  before.  Ihere  Is  evid¬ 
ence  bovever,  ^Ich  Indicates  that  particles  are  continuing  to  crush 
throughout  Region  4.  Chapter  4  considers  the  energy  required  to  create 
new  surfaces.  It  Is  found  that  after  crushing  begins  the  energy  absorbed 
due  to  particle  crushing  seems  txj  be  linearly  related  to, the  stress  level. 

Unfortunately,  there  are  no  data  available  in  the  literature 
which  can  be  used  to  fully  corroborate  the  analytical  conclusions  sum¬ 
marized  above  and  presented  In  detail  in  Chapters  3  and  4.  In  order  to 
fill  this  gap  In  our  experimental  knowledge  a  series  of  exploratory  tests 
on  smooth  uniform  sand  In  one  dimensional  compression  were  conducted  to 
provide  a  better  understanding  of  the  mechanisms  by  which  sand  absorbs 
energy.  The  design  and  construction  of  a  new  testing  apparatus  was  re¬ 
quired  to  obtain  the  necessary  data  In  the  higher  pressure  ranges.  The 
design  of  the  apparatus  and  the  test  results  are  described  In  Chapter 

The  data  presented  in  Chapter  5  show  that  energy  can  be  absorbed 
by  sand  In  stress  ranges  below  crushing  with  only  very  small  permanent 
strains  resulting.  This  behavior  is  substantially  the  same  as  predicted 
from  the  analysis  of  the  rnedlun  composed  of  equl -radii  spheres. 
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CHAPTER  2 

SEIfCTED  SUMAFX  OF  PREVIOUS  WORK 
It  is  not  the  puz^ose  herein  to  make  a  thorough  survey  of  all 
of  the  Important  restilts  reported  In  the  literature  vblch  have  contributed 
to  an  understanding  of  the  behavior  of  granulcu:  materials  and  soils.  Such 
an  effort  wotild  be  quite  voluminous  and  an  extensive  report  in  Itself.  It 
Is  the  writers'  tmderstandlng  that  Professor  R.  V.  Whitman  at  MIT  Is  cur¬ 
rently  preparing  such  a  report,  [2.371*  Furthermore,  a  very  excellent 
comprehensive  svimmary  of  the  state  of  knowledge,  on  the  behavior  of  granu¬ 
lar  materials,  up  to  1937  has  been  made  by  Deresiewicz  [2.ll.  !nie  reader 
is  referred  to  these  works  for  a  comprehensive  survey  of  the  literature  on 
the  behavior  of  soils. 

The  purpose  of  this  chapter  is  to  call  attention  to  some  of  the 
more  recent  findings  uhlch  are  important  In  understanding  soil  behavior 
and  also  to  point  out  some  other  works  irtilch,  while  not  new,  have  been 
someidiat  overlooked  In  the  literature.  It  Is  hoped  that  the  few  comsents 
made  here  will  be  useful  to  the  reader  as  supplements  to  the  two  above 
noted  works. 

2.1  Theories  of  Granular  Media 

There  are  many  theories  which  have  been  developed  In  the  past 
to  attempt  to  predict  or  rationalize  the  behavior  of  granular  materials. 

Some  of  these  have  been  based  on  the  concepts  of  the  theory  of  plasti¬ 
city.  Others  have  been  foxmded  on  a  study  of  spherical  arrays  such  as 
that  sunarlzed  by  Deresiewicz  [2.1]. 

A  recent  publication  worthy  of  note  Is  a  susaary  of  Russian 
contributions  In  soil  mechanics  by  Klein  [2.2].  This  Is  primarily  a 
historical  account  of  the  developswnt  of  the  theory  of  granular  media  In 

I 


Russia  from  1917  to  1957  aa<i  includes  little  actual  theory.  Its  major 
contribution  Is  that  It  provides  an  excellent  bibliography  of  only  Rus¬ 
sian  vork  (some  175  listings)  in  this  area  up  to  the  early  part  of  1957* 
Another  survey  of  recent  Russian  work  In  soil  mechanics  was  also  con¬ 
ducted  by  Drashevska  la.jVj  in  1958* 

One  of  the  most  recent  theories  of  soil  mechanics  developed 
In  Russia  Is  that  presented  by  Grigoryan  [2.5]  >  He  formulates  a  the¬ 
ory  ^Ich  Is  presumed  to  be  applicable  for  soils  subjected  to  extremely 
large  pressures  (in  terms  of  thousands  of  psl).  His  theory  suggests 
that  the  behavior  of  a  soil  Is  dependent  on  two  functions^  one  character¬ 
izing  the  volvsne  deformation  of  the  medium  and  the  other  the  properties 
of  the  medium  during  shear  deformation.  It  should  be  noted  that  the 
most  important  difference  between  the  Grigoryan  model  and  the  usual 
plasticity  theory  is  that  here  the  volume  deformation  relationships  be¬ 
tween  the  Man  pressure  p  •  -l/?  (o  +  a  +  c  )  and  the  density  la  not 

X  y  z 

assumed  to  be  elastic  (reversible).*  Ihey  are,  however,  connected  by  a 
well  defined  relationship  idilch  Is  different  for  Increasing  and  decreas¬ 
ing  pressure.  This  theory  gives  a  loading  and  unloading  stress-strain 
curve  similar  In  shape  to  Fig.  2.1.  Grigoryan's  theory  Is  discussed  In 
greater  detail  In  a  later  publication  [2.k]  and  bom  Inaccuracies  In 
the  original  work  eire  corrected. 

SoM  experlMntal  work  Is  presented  by  Alekseenko,  Grigoryan, 
et  al  [2.4]  In  support  of  the  theory  proposed  by  Grigoryan  for  granulcu: 
media.  A  series  of  TVT  explosions  was  set  off,  and  the  soil  behavior 
ma  Masured  radially  outward  away  from  the  blasts.  The  results  seeMd 
to  show  a  definite  linear  relationship  between  the  Man  compressive  stress 
and  the  square  root  of  the  second  stress  Invariant.  For  average  stresses 
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les8  than  15  kg/ cm^  It  was  found  that  the  second  stress  Invariant  could 
be  related  to  the  average  stress. 

Another  recent  theory  of  soil  mechanics  has  been  presented  In 
a  monograph  by  Geniev  [2.5I.  Here  a  theory  for  the  plane  motion  of  a 
granular  medium  Is  formulated  \dilch  is  an  extension  of  the  work  of  Soko¬ 
lovsky.  The  equations  of  motion  consist  of  (l)  two  equations  of  equl- 
llbrlm  (2)  conditions  for  Coulomb  boundary  equilibrium  (3)  Incompress¬ 
ibility  conditions  and  (^)  conditions  of  convergence  of  the  directions 
of  maximum  velocities  of  shear  deformation  vlth  the  direction  of  lines 
of  slip.  The  application  of  Geniev's  theory  to  the  dynamic  soil-struc¬ 
ture  Interaction  problem  Is  currently  under  study  at  the  University  of 
Illinois  [2.35]. 

One  theoretical  treatment  ^Ich  is  not  too  veil  knovn  is  that 
developed  by  Oshima  [2.6]  vherein  an  attempt  Is  made  to  specify  in  tensor 
notation  all  of  the  contributing  effects  vhlch  should  be  incltided  in  a 
theory  of  gran\ilar  media.  The  approach  has  some  merit  In  that  an  attempt 
Is  made  to  present  a  consistent  theory  which  encompasses  the  entire  geo¬ 
metry  of  the  medium.  Ghfortunately,  while  the  abstractness  presents  a 
convenient  compact  form  for  the  theory,  there  Is  difficulty  in  obtaining 
all  of  the  terms  In  the  significant  tensors. 

Katz  and  others  [2.7  aad  2.8]  have  studied  the  propagation  of 
plane  stress-waves  In  sand  with  a  two-parameter  exponential  stress-strain 
lav,  Including  hysteresis.  Along  a  similar  line  Parkin  has  treated  the 
one  dimensional  wave  propagation  In  sand  in  an  atteiqpt  to  explain  the  be¬ 
havior  measured  by  Whitman  [2.I3]  in  some  experiments  carried  out  In  the 
early  1930' s.  The  subsequent  discussions  of  the  latter  paper  by  Fulton 
and  Hendron  [2.II],  Sellg  and  Vey  [2.12],  Whitman  [2.I3]  and  Parkin  [2.14] 
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are  particularly  interesting  in  euBnarlzlng  the  current  state  of  know¬ 
ledge  and  problems  associated  with  predicting  the  behavior  of  sand.  An 
investigation  of  the  dynamic  soil -structure  interaction  problem  using  a 
discrete  mass -spring  model  to  represent  the  soil  is  currently  under  way 
at  the  University  of  Illinois  I2.56]. 

Some  of  the  more  recent  papers  in  English  on  the  application 
of  the  theory  of  plasticity  to  the  behavior  of  soils  include  [2.I6  - 
2.20].  Hieoretlcal  studies  of  the  stress  conditions  existing  in  tri- 
axlal  compression  have  been  formulated  by  Balia  [2.2l]  and  Haythomthwalte 

[2.22]. 

2. 2  Experimental  Results 

One  result  ^ich  is  of  Importance  to  this  study  is  given  in  a 
paper  by  Allen,  et  al  [2.25],  vdilch  contains  the  results  of  a  series  of 
tests  on  sand  (Fig.  2.2).  The  particular  significance  of  these  tests  is 
that  the  pressures  were  extremely  hig^i  (approximately  90,000  psl).  Ilhe 
sample  was  placed  inside  a  steel  cylinder  to  approximate  conditions  of 
one  dimensional  compression.  Unfortunately,  the  sample  size  was  of  the 
wrong  proportion  for  a  one  dimensional  test  (2.31  inches  in  diameter  and 
7.66  Inches  high)  BQd  the  stress-strain  behavior  may  very  well  be  signi¬ 
ficantly  affected  by  the  deformations  of  the  container  and  side  friction 
effects.  Ihese  results  are  of  significance,  however,  in  indicating  the 
general  shape  of  the  stress -strain  curve  in  the  high  pressure  regions. 
Grigoryan  used  this  general  shape  to  support  his  theory  of  granular  media. 

J.  M  Roberts  [2.24]  has  investigated  the  hysteresis  character¬ 
istics  of  an  Ottawa  sand  with  lateral  confinement  closely  approximating 
one  dimensional  conditions.  Some  200  loading  cycles  were  carried  out  in 
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a  stress  region  belov  50  psi.  It  was  found  from  the  hysteresis  loops  that 
sand  can  absorb  energy  without  having  any  permanent  displacements.  The 
reason  for  this  phenomenon  noted  by  Roberts  and  more  conclusive  experi¬ 
mental  results,  are  given  by  the  authors  In  the  present  report. 

Studies  on  the  one  dimensional  stress-strain  behavior  of  sand 
and  quartz  in  the  higher  pressure  regions  were  carried  out  by  J.  E.  Rob¬ 
erts  and  others  at  MIT  [2.25,  2.26].  Some  of  these  data  are  analyzed  in 
Chapter  U  of  this  report  in  an  effort  to  shed  some  light  on  the  mechanism 
by  idilcb  sand  absorbs  energy  in  pressure  regions  of  ^,000  -  20,000  psi. 

Tests  on  the  static  and  dynamic  stress -strain  characteristics 
of  sands  in  the  low  pressure  regions  have  been  conducted  by  Thompson  at 
the  Ballistics  Research  Laboratory  and  the  results  have  been  published 
in  a  series  of  short  reports  [2.27,  2.28,  2.29J. 

KJellman  and  JWcobson  [2.50]  investigated  the  influence  of  grain 
shape  and  grain  size  on  the  behavior  of  granular  materials  incased  in  a 
series  of  steel  rings  in  an  attempt  to  eliminate  the  side  friction  effects. 
A  trlaxial  compression  apparatus  was  modified  by  Welssman  and  Hart  [2.31I 
to  study  the  damping  capacity  of  granular  soils  as  reflected  by  the  hys¬ 
teresis  loop.  Wolfskin  and  Buchanan  [2.32]  discuss  some  of  the  aspects 
of  dynamic  foundation  loads  supported  by  a  granular  soil.  Other  recent 
work  on  dynamic  trlaxial  tests  includes  that  by  Shannon,  et  al  [2.33]* 
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p^=  Initial  Mean  Stress 
dp-  Initial  Density 


Fig.  2.  1  ASSUMED  STRESS  STRAIN  CURVE  FOR  SAND 
(after  Grigoryan  (2.  3,  2.4)) 


Strain 

Fig.  2.2  ACTUAL  STRESS  STRAIN  CURVE  FOR  SAND 
(after  Allen,  et  al  (2.  23)) 
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CHAmR  5 

THSORETICAL  APPROACH  FOR  PREDICTDIG  THB  STRESS-STRAIN 

BEHAVIOR  OF  ISNSE  SAHBS  PRIOR  TO  PARTICLE  CRUSHING 

The  stress-strain  behavior  of  most  materials  In  the  stress 
levels  of  Interest  may  be  assumed  to  approximately  follow  Hooke's  lav.  If 
the  materlaO.  Is  isotropic,  the  stress-strain  relations  can  be  completely 
described  by  tvo  independent  elastic  constants. 

nie  discussion  in  the  introduction  to  this  report  pointed  out 
the  fact  that  constants  could  not  be  used  to  describe  the  stress -strain 
properties  of  a  granular  medium.  It  vas  shown  that  it  is  possible  to 
obtain  almost  any  stress-strain  curve  one  desires  by  varying  the  stress 
level,  the  state  of  stress,  and  displacement  conditions  at  the  boundaries. 
Thus,  there  is  no  unique  stress -strain  relation  for  a  granular  siedium. 

This  does  not  mean,  however,  that  the  stress -strain  relations  of  certain 
prescribed  granular  media  under  a  given  set  of  boundary  conditions  can¬ 
not  be  found.  The  above  merely  indicates  that,  even  for  a  greatly  ide^- 
ized  medium  such  as  a  uniform  array  of  equl-radll  si^eres,  the  problem 
of  defining  a  stress -strain  relation  for  a  granular  medium  is  a  boundary 
valued  problem  in  particulate  mechanics  and  any  euialytlcal  study  of  be¬ 
havior  under  changes  in  stress  must  adequately  take  into  account  the  bound¬ 
ary  conditions. 

A  Biecbanics  approach  for  formulating  the  behavior  of  this  type 
of  media  must  consider: 

1.  The  equilibrium  of  each  particle  and  the  medium  as  a  ^ole. 

2.  Certain  geometrically  admissible  conditions  on  the  deforma¬ 
tions  of  the  particles  (i.e.,  a  set  of  compatibility  equa¬ 
tions). 

3.  The  relationships  between  the  normal  forces  and  normal  dis¬ 
placements  and  shearing  forces  and  tangential  displacements 
at  each  contact  point  on  a  particle. 
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4.  Hie  boondary  eoadltlone  on  the  MdlioD. 

Ihe  non-linear  partial  differential  equations  idilch  result 
froBi  even  the  slnplest  aodel  of  equl-radll  spheres  are  extreaely  complex 
and  closed  form  solutions  can  be  obtained  for  only  a  fev  cases.  Since 
some  of  the  solvable  problems,  hovever,  turn  out  to  be  basic  ones,  such 
an  approach  has  a  great  deal  of  merit  In  adding  to  the  understanding  of 
the  behavior  of  granular  materials. 

This  chapter  revlevs  the  formulation  of  the  theory  of  Granular 
Hedla  eccording  to  the  ideas  of  Mlndlln  and  Suffy  and  then  presents  a  modi¬ 
fication  of  the  theory  vhlch  is  appropriate  for  one  dimensional  compress¬ 
ion. 

3,1  The  Mlndlln -Duffy  Theory  of  Granulcu-  Media 
3.1.1  Basic  Hieory 

The  theory  developed  here  follows  that  first  formulated  by 
Duffy,  Mlndlln  (3-1,  3.2,  3.3J  and  their  co -workers.  Hie  granular  med¬ 
ium  is  restricted  to  a  face  centered  array  of  equl-radll  spheres  (Fig. 

3.1).  The  face  centered  array  was  chosen  primarily  because  this  eurray 
is  one  vhlch  provides  for  the  densest  possible  packing.  Hence,  there  Is 
no  vay  In  vhlch  compressible  deformations  of  the  medium  can  take  place 
as  a  result  of  rearrangement  of  the  particles. 

Consider  the  behavior  of  two  spheres  pressed  together  by  a 
normeO.  force.  According  to  the  Hertz  theory,  if  two  elastic  spheres, 
each  having  radius  R,  shear  modulus  M,  and  Poissons  ratio  v,  are  mut¬ 
ually  compressed  by  a  normal  force  N,  the  resulting  surface  of  contact 
Is  a  plane  bounded  by  a  circle  of  radius 

1/3 
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and  the  rate  of  change  of  the  relative  approach  of  the  sphere  centers, 
a,  vlth  respect  to  the  normal  force  If,  Is 


C  =  ^ 

dU  2Ma 


(5.2) 


where  C  Is  defined  as  the  normal  compliance.  Two  spheres  In  normal  con¬ 
tact  are  shown  in  ?lg.  3.2  while  Fig.  3.3  shows  the  distribution  of  the 
normal  stress  across  the  area  of  contact.  A  complete  derivation  of  the 
Hertz  compliance  and  the  assumptions  and  limitations  of  the  Hertz  theory 
In  regard  to  this  study  appear  In  Appendix  A. 

Suppose  now  that  the  system  of  like  spheres,  Initially  compressed 
by  a  constant  normal  force  H  Is  subjected  to  a  tangential  force  T  which 
acts  In  the  plane  of  contact  and  whose  magnitude  Increase  monotonlcally 
from  zero  to  a  given  value.  Because  of  symmetry,  the  distribution  of 
normal  pressure  remains  imchanged.  If  It  is  assumed  that  there  Is  no 
slip*  on  the  contact,  then,  because  of  symmetry,  the  displacement  of  the 
contact  sxirface  in  Its  plane  Is  that  of  a  rigid  body.  Die  solution  of 
the  appropriate  boundary -value  problem,  due  to  Cattaneo  [3.4)  and  Mlnd- 
lln  [3.5J>  yields  the  tangential  component  of  traction  t  on  the  contact 
surface  and  the  tangential  displacement,  &,  of  points  In  one  sphere  re¬ 
mote  from  the  contact  vlth  respect  to  similarly  situated  points  In  the 
other  sphere.  Die  tangential  traction  Is  parallel  to  the  applied  force 
T,  axially  symmetric  In  magnitude,  and  increases  without  limit  on  the 
bounding  curve  of  the  contact  area.  Fig.  3.3.  It  is  reasonable  to  sup¬ 
pose  that  slip  Is  initiated  at  the  edge  of  the  contact  since  It  Is 
there  the  singularity  In  traction  takes  place  in  the  absence  of  slip. 

Since  without  slip  this  traction  Is  symmetric,  the  slip  Is  assumed  to 


*  By  slip  we  mean  relative  displacement  of  contiguous  points  on  a  portion 
of  the  contact  surface.  We  distinguish  between  "slip”  axid  the  term 
"sliding,"  which  we  reserve  to  denote  relative  displacement  over  the  en¬ 
tire  contact. 
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progress  radially  Inward,  covering  an  annular  area.  On  this  annulus  It 
Is  assumed,  as  a  first  approximation,  that  the  tangential  component  of 
traction  Is  In  the  direction  of  the  applied  force  and  Is  related  to  the 
normal  (Hertz)  component  of  stress  a,  already  present.  In  accordance 
with  Coulomb's  Lav  of  sliding  friction.  Hence: 

T  -  fo  (3.5) 


where  f  Is  the  coefficient  of  static  friction  and  a  Is  given  by 


a 


(3.3a) 


Here  p  Is  the  distance  .from  the  center  of  the  contact  circle. 
Hie  distribution  of  t  with  slip  is  also  shown  in  Fig.  3 •5.  Cattaneo  and 
Mlndlln  also  predict  the  relation  between  the  radius  of  the  adhered  por¬ 
tion  and  the  applied  tangential  force  as 


b 


1/3 


(3.4) 


where  b  is  as  shown  in  Fig.  3*3  and  where  the  tangential  compliance  at 
the  contact  Is  given  by 


S  = 


12. 

dT 


2  -  V 

W 


-l/3 


(3.5) 


Equation  (3.3)  bolds  only  for  the  case  of  an  Increasing  T  with 
a  constant  N  at  the  contact.  In  the  problems  we  will  consider  in  this 
paper,  T  and  H  are  both  Increasing  so  the  tangential  compliances  appropri¬ 
ate  for  this  problem  are  given  In  reference  [3.6J  as 
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S  = 


2  -  V 

w 


(3.6) 


or 


S  = 


2  -  V 

W  ' 


Vf 

dT  - 


(3.7) 


Let  us  Dov  coQsider  a  face  centered  cubic  array  of  uniform 
spheres  as  shown  in  Fig.  3>1>  ^ere  each  sphere  Is  In  contact  with  12 
other  spheres.  A  typical  element  of  this  packing  Is  shown  in  Fig.  3.U. 
If  the  coordinate  system  shown  In  Fig.  3*^  Is  translated  to  the  center 
of  any  sphere.  Fig.  3-5>  the  12  contact  points  woxxld  have  the  following 
coordinates  in  terms  of  the  radius  of  the  sphere. 
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There  are  two 

shearing 

forces  and 

one  normal  force 

at  each  con 

tact  point  as 

shown  on  Fig.  3.5. 

In  Index 

notation 

It  Is  convenient  to 

Identify  the  various  components  of  normal  forces  by  the  symbol. 


\rtiere  the  subscripts  correspond  to  the  planes  In  which  the  components  lie. 


The  two  tangential  force  components  at  each  contact  are  chosen  to  lie  In 
and  normal  to  the  co-ordinate  pleuses,  and  are  Identified  by  and  ’itk 
respectively.  At  all  contacts  ^ere  the  normal  has  direction  cosines  of 


unlike  sign  the  force  components  are  further  distinguished  by  primes. 
Thus  the  components  at  the  contact  points  1  through  12  are : 
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Let  UB  nov  consider  the  equlllbrlvoD  equations  of  a  representa¬ 
tive  cube  (Fig.  3*^)  due  to  an  arbitrary  Initial  state  of  stress  imposed 
on  the  medium.  If  an  Increment  in  stress  is  added  to  the  initial  stress, 
the  forces  on  the  faces  of  the  cube  vlll  be  increased  by  corresponding 

increments  dP  .  Let  us  determine  the  unknown  increments  of  force  at  the 
^  J 

contacts  between  the  spheres  within  the  cube  idiich  result  troa  the  force 
increments  dP  .  If  the  stress  increment  in  the  medium  is  homogeneous, 
then  the  contact  forces  will  be  equal  at  contacts  having  corresponding 
positions  on  the  surfaces  of  the  sj^eres.  Ibe  contact  forces  diametri¬ 
cally  opposed  on  each  sphere  are  equal,  thus  only  l6  of  contact  forces 
on  each  sphere  are  independent.  Hot  only  must  the  cube  as  a  whole  remain 

in  equilibrium  under  the  action  of  the  Increments  in  applied  force,  dP  , 

^  J 

but  each  sphere  and  each  portion  of  a  sphere  within  the  cube  must  also 
remain  in  equilibrium.  Since  the  portions  of  spheres  in  the  cube  axe 
acted  on  both  by  applied  forces  and  by  contact  forces,  the  equations  of 
equilibrium  relate  the  Incrssents  in  applied  forces  to  the  increments  in 
contact  forces.  By  writing  equilibrium  equations  for  various  octants  of 
spheres  in  the  cube  as  shown  in  Fig.  3*6  it  may  be  shown  that  there  are 


nine  Independent  equations  of  equilibrium,  niese  equations  are  the  fol¬ 


lowing: 
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The  nine  equilibrium  equations  given  above  are  not  sufficient 
to  uniquely  describe  the  behavior  of  the  medium  since  the  problem  is 
statically  indeterminate.  A  set  of  compatibility  conditions  are  formu* 
lated  by  giving  consideration  to  the  admissible  displacements  of  the 
medium  vhich  result  frcxn  the  incremental  stresses.  Let  the  components 
of  relative  displacements  of  the  centers  of  spheres  be  designated  by 
dQ^j,  ^  correspond,  respectively,  to  the  forces 

dT  ,  dT,,  (Fig.  3*7).  We  will  require  that  the  displacements  of  the 

X  j 

center  of  the  spheres  be  single  valued,  l.e.,  the  vector  distance  auK>und 
the  closed  path  through  the  center  of  the  spheres  must  vanish  both  be¬ 
fore  and  after  the  medium  is  strained.  Hence  the  sum  of  the  relative 
displacements  of  the  center  of  these  s];heres  around  the  closed  path  must 
vanish.  Expressions  for  this  condition  for  all  possible  paths  connect¬ 
ing  the  center  of  a  sphere  in  the  medium  yields  9  Independent  equatloxis 
of  compatibility  as  follova: 
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(3.9) 


Ibe  compatibility  equations  can  be  vrltten  In  terms  of  force 
Increments  by  using  the  compliances 


da'  -  C  dJT' 

"ij  ■ 

(3.10) 

"V  ■  WW 

Where  the  subscripts  Identify  the  contact.  Ibus  Eq.  (5*9)  together  vlth 
Eq.  (3.3)  yield  18  independent  equations  containing  the  applied  forces 
dPj^j,  and  l8  independent  ccmiponents  of  the  unknovn  contact  forces, 
dTj^j,  and  dT^.  Unfortunately,  hoverer,  these  18  equations  are  non-lin¬ 
ear  since  the  equations  Include  the  compliances  and  the  compliances  are 
themselves  functions  of  the  contact  forces.  Consequently,  a  solution  to 
these  equations  Is  a  very  difficult  problem  and  can  be  obtained  only  for 
certain  simple  cases. 

Ihe  incremental  extensional  and  shearing  strains  in  the  array, 
expressed  in  terms  of  the  compliances,  are 


(3.11) 


Llkevlse  the  applied  force  Increments  in  Fig.  are  related  to  the 
stress  Increments  for  a  face  centered  array  by 
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dPjj  .  (5.12) 

3.1.2  Application  to  a  gydrostatlc  State  of  Strees 

One  problem  of  Interest  vhlch  can  be  solved  Is  that  of  a  face 
centered  array  subjected  to  a  hydrostatic  state  of  stress.  Under  this 
state  of  stress  we  can  obtain  from  symmetry  the  conditions  that  all  norm¬ 
al  forces  at  the  contacts  are  equal  and  all  shearing  forces  are  equal. 

If  we  let 

P  =  total  force  on  the  face  of  a  differential  element 

o 

N  «  normal  contact  force 
o 

=  tangential  contact  force 
R  =  radii  of  spheres 
the  equilibrium  Eqs.  (3*8)  reduce  to 

4dT^  +  W2dJI^  =  dP^  (3.15) 

4dT  -  W2d»  ■=  dP 
o  00 

Adding  Eqs.  (3* 13)  and  considering  the  InltleQ  condition  of  zero  stress 
yields  the  expected  condition  for  the  tangential  force 

T  »  dT  »  0 
o  o 

Hence,  the  equilibrium  equation  becomes 


<^0 


If  Is  the  hydrostatic  strain  in  any  direction,  Eq.  (3 *11) 

yields 


(3.15) 


-27- 


vhere  da^  Is  the  normal  deformation  due  to  dJI  and  is  related  to  dM  by 
o  o  0 

the  normal  compliance  C^. 


da  «  C  dl 
o  0  o 


(3.16) 


From  the  Herts  theory 


C 

o 


1  -  V 

2Ma^ 

o 


(3.17) 


vhere  a  la  the  radius  of  contact  defined  by 


p3(i  -  y)  V  1 

‘L - - J 


V3 


(3.18) 


By  combining  Sqs.  (3.I5),  (3.16),  (3.17)  and  (3.18),  the  Incremental 
strain  may  be  expressed  as 


1/2 


r  r-.2/5 


(3.19) 


Integrating  Eq.  (3.19)  and  considering  zero  Initial  conditions  results 


In 


'3(l-v)  V2' 

8I» 

W  a 


?/3 


(3.20) 


It  should  be  noted  that  Eq.  (3.20)  Is  also  the  equivalent  of  the  first 
stress  Invariant  I,  being  proportional  to  the  three  halves  pover  of  the 
first  strain  Invariant  J. 

It  Is  of  Interest  to  cotspare  the  results  of  this  theory  with 
some  limited  experimental  results  available  In  the  literature.  KJellmax- 
[3.7]  carried  out  some  tests  on  dry  sand  subjected  to  hydrostatic  pres¬ 
sure  and  measmred  the  strains  associated  with  the  pressure.  Table  3.I 
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gives  the  variation  of  the  first  stress  and  strain  invariants  taken  from 
his  experimental  results.  A  plot  of  this  Is  shorn  in  Fig.  ^.6,  Also 
shovn  is  a  plot  of  the  exponential  relationship  indicated  by  the  theory 
discussed^  vhlch  is  fitted  to  the  experimental  curve  at  the  6  kg/cm 
stress  level.  Ibe  behavior  of  the  sand  is  not  as  stiff  as  the  dense 
packing  theory  Indicates  in  the  lover  pressure  regions  while  it  becomes 
stiffer  than  the  theory  predicts  in  the  higher  pressure  regions. 

The  former  seems  reasonable  in  that  in  the  lover  pressure  reg¬ 
ion  the  number  of  contacts  is  not  as  large  as  the  12  assinaed  for  a  dense 
state.  Hence  the  real  sand  is  not  as  stiff  as  the  spheres.  In  the  hl^- 
er  pressure  regions ,  it  seems  reasonable  to  associate  the  stiffening 
effect  vlth  the  random  size  distribution  of  particles  vhlch  differs  from 
the  equl-radll  spheres  assumed  in  the  theory.  If  the  sand  grains  are 
of  different  sizes,  more  than  12  contacts  can  be  reali.zed  after  the  sand 
has  been  compressed  sufficiently  to  allow  these  contacts  to  be  made. 

This  would  tend  to  make  the  sand  stiffer  than  the  spheres.  Furthermore, 
while  every  attempt  was  made  by  KJellman  to  minimize  the  effect,  any 
frictional  resistance  due  to  the  testing  apparatus  would  effectively 
stiffen  the  measured  stress -strain  relations. 

From  the  above  discussion  and  development,  it  can  be  concluded 
that  the  theory  associated  with  the  face  centered  array  of  spheres  shows 
at  least  some  qualitative  correlation  with  the  behavior  of  dense  sand 
subjected  to  a  hydrostatic  state  of  stress. 

5.2  One  Dimensional  Theory  of  Granular  Media 
5.2.1  Monotonlcally  Increeislng  Load 

Let  us  now  extend  the  theory  of  Dufiy  and  Mlndlln  developed  in 
Section  5.1.1  to  the  solution  of  the  stress -strain  behavior  of  eui  array 
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of  spheres  subjected  to  one  dimensional  compression.  One  dimensional 
compression  is  defined  as  that  state  of  stress  resulting  from  the  appli¬ 
cation  of  a  load  in  the  vertical  or  2  direction  when  the  lateral  strains 
in  the  x  and  y  directions  are  zero. 

The  radial  symmetry  of  the  one  dimensional  problem  greatly 
simplifies  Eqs.  (3-8)  and  (5-9)-  The  representative  cube  now  beconies 
Fig.  3-9>  and  the  forces  on  a  sphere  reduce  to  those  shown  on  Fig.  3. It. 
From  symmetry  the  following  simplifications  can  be  made  for  the  forces 
and  displacements ; 
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the  associated  compliances  now  become 
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Using  the  above  simplifications  and  considering  symmetry  the 
equilibrium  Eq.  (3.8)  reduce  to 


\l2 

dUg  +  dTg  =  dP 


zz 


7  2 

dN^  +  dNg  -  dTg  =  dPg 


(3.24) 


In  a  similar  fashion  the  compatibility  Eq.  (3.9)  reduces  to 


dOi  -  da^  +  dbg  =  0 


(3.25) 


wner 


e  rdi'j  compliance  Eqs.  (3.IO)  now  become 


dai  =  C^dN^ 


da  =  C  dN 
2  2  2 


(3.26) 


dOg  =  SgdTg 


Substituting  Eqs.  (3*26)  into  Eq.  (3.25)  yields 


CidNi  -  CgdNg  +  SgdTg  =  0 


(3.27) 


Equations  (3.24)  and  (5.27)  axe  sufficient  to  describe  the  be¬ 
havior  of  a  granular  medium  subjected  to  given  vertical  and  lateral  forces 

P  and  r„.  If,  however,  only  the  vertical  force  is  known  and  the  lateral 

Z  u  n 

force  must,  also  be  determined,  a  further  condition  is  necessary.  This 
condition  naturally  comes  from  the  lateral  strain  relationship. 

Due  to  symmetry  the  lateral  strain  c„  determined  from  Eqs.  (5.II) 

li 


reduces  to 


(5.28) 
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or 


‘^"h  = 


2R 


(3.29) 


Thus  the  behavior  of  a  face  centered  array  of  spheres  subjected 
to  a  vertical  force  and  restricted  to  symmetrical  lateral  deformations 
can  be  obtained  from  a  solution  to  the  following  equations: 
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(3.30) 


2Rd£ 


dlf. 
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For  the  case  of  Interest  here,  namely  one  dimensional  compress¬ 
ion,  ve  require  €„  to  vanish  for  all  loadings.  For  this  particular  stress 
a 

state  euid  zero  initial  conditions,  we  obtain  =>  0  and  Eqs.  (3 '30)  re¬ 
duce  to 


/ 2 

dRg  +  (iTg  .  -g-  dP^^ 

(3.31a) 

sj  2 

(iNg  -  dTp  -  ip  dPg  =  0 

(3.31b) 

s 

dNg  -  dT^  =  0 

(3.31c) 

2 

^ere  the  expression  for  —  is  obtained  from  Eqs.  (3-6)  and  (J.T)  as: 

^2 


r  dNp  dNp  T  -1/3-1 

■  ‘i  5^  *  f  1  -  ?»;>  J 


where 
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2-v 


(5.32b) 


Hie  compatibility  Eq.  (3.31c)  states  that 


dOg  =  dB^ 


Hence  for  the  one  dimensional  case  Eq.  (3.32)  reduces  to 


da  c  dfl 
2  2  2 


"zz  R 


(3.33) 


It  should  be  noted  that  Eqs.  (3.3l)  are  coupled  non-linear 

differential  equations  because  of  the  compliances  in  the  third  equation. 

The  vertical  strains  e  associated  with  the  behavior  of  this 

zz 

medium  can  be  obtained  from  Eqs.  (3.1l)  as 


(3.3^) 


For  the  solution  of  Eqs.  (3.3l)i  let  us  consider  first  the  com¬ 
patibility  Eq.  (3.31c).  Substituting  Eqs.  (5.52)  into  Eq.  (3.31c)  yields 
after  some  rearrangement  of  terras 


1-k  f  T 

f  >  (1  -  ^) 

Ki  rrtg 


(3.35) 


Introduce  now  the  new  veu"iable 


(3.36a) 
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NgZ^  =  Ngf  -  Tg 

dNgZ^  +  )  Z^dZ  =  fdNg  -  dTg 

or 

dT 

55^.  f  -  z5  ,  .  (3.56b) 

Using  Eqs.  (3.36),  Eq.  (3.35)  can  be  transformed  into  the  fol¬ 
io  ving  form. 


^ere 


3zaz 

2 

Z  +K^ 


1-k^f 


Integrating  both  sides  results  in 


In  1*2  -  In  A  -  -3/2  In  (Z^  +  Kg) 


or 


T-  <^‘^^"2) 


-3/2 


vhere  A  is  a  constant  of  integration. 
From  Eq.  (3»36a),  Eq.  (3.39)  becomes 


(3.37) 


(3.38) 


(3.39) 


(3.40) 


Equation  (3-40)  is  the  general  solution  to  Eq.  (3.31c)  and  a 
particular  solution  may  be  obtained  by  evaltiating  the  constant  of  integration 
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A.  Since  we  are  Interested  in  the  initial  conditions  P  =  =  0, 

one  is  Inclined  to  evaluate  A  from  these  conditions.  While  this  is  a 
true  boundary  condition^  a  singularity  point  occurs  at  the  origin.  Hie 
equation  is,  however,  well  behaved  at  other  points  and  the  following 
tecnnique  will  be  used  to  evaluate  A. 

Imagine  a  small  hydrostatic  state  of  stress  initially  holding 
the  spheres  in  contact  before  the  one  dimensional  state  is  Imposed.  This 
hydrostatic  stress  produced  no  initial  tangential  force  at  the  contacts 

but  it  does  cause  a  normal  force  of  =  N  . 

2  0 

From  Eq,  (j. 40) 

N  =  A(f2/3  + 
o  2 


or  solving  for  A,  the  constant  of  integration 

N 


Substituting  A  into  Eq.  (3 
and  taking  into  account  Eq,  (3.38)> 


(3.41) 


.40)  yields  after  some  rearranging 


T 


2 


N  2/3 


(— ) 


(3.42) 


A  plot  of  a  family  of  curves  representing  Eq.  (3*42)  with  var¬ 
ious  values  of  the  initial  hydrostatic  stress  is  given  in  Fig.  3.11. 
Hie  paths  of  loading  are  indicated  by  arrows  on  the  curves.  For  conven¬ 
ience  the  coefficient  of  friction  f  was  taken  as  O.3  and  Poisson's  ratio 


V  as  0.2. 
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If  an  initial  hydrostatic  stress  is  imposed  we  can  determine 
the  value  of  when  the  sliding  occurs.  For  this  condition  we  see  that 
from  Eq.  (3.42) 

r  “n  2/5  -I 

[(^)  -  1  ]  .  1  =  0 

Rewriting  euid  taking  account  of  (3.32h)  yields 


(3.45) 


The  reader  will  recall  that  we  were  Initially  interested  in 
defining  the  relation  between  and  for  a  granular  medium  in  a  one 
dimensional  state  of  stress  which  was  loaded  from  an  unstressed  condi¬ 
tion.  We  noted  that  the  solution  contained  a  singularity  at  the  zero 
stress  and  it  was  necessary  to  provide  an  initial  hydrostatic  stress 
to  hold  the  spheres  in  contact  when  loading.  If  we  now  allow  to 
approach  zero  we  obtain  the  solution  desired. 

On  investigating  Fig.  3-11  or  Eq.  (3.42)  it  is  clear  that  if 
one  dimensional  loading  commences  from  a  completely  unstressed  condition 
(n^  =  0),  the  relation  between  and  Tg  is  that  of  a  straight  line  with 

a  slope  of  l/f.  Thus,  the  relationship  between  Ng  and  Tg  for  a  face 
centered  array  subjected  to  one  dimensional  compression  is 

T  =  fH 

2  2 

This  means  that  sliding  (differentiated  from  slip  in  Section 
3.1.1)  at  the  contacts  occurs  immediately  on  initiation  of  loading.  This 
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is  extremely  important  in  that  the  theory  developed  to  date  along  the 
lines  of  Mindlin>  et  al,  does  not  allow  sliding  at  the  contacts.  Ihe 
authors^  however,  will  show  in  the  following  that  the  stress-strain  rela 
tion  ships  in  the  one  dimensional  granular  medium  ccm  he  described  mathe 
matically  even  though  sliding  occiirs  throughout  the  application  of  a 
monotonically  increasing  load. 

From  compatibility,  even  in  the  case  of  sliding,  the  geometri¬ 
cal  relationship  still  holds  for  the  displacements  at  the  contacts,  l.e. 

dOg  -  dbg  -  0  (5.44) 

It  should  be  clear,  however,  that  the  tangential  displacement 
6^  is  no  longer  related  to  the  tangential  contact  force  Tg  by  the  tangen 
tlal  compliance  Sg  because  of  sliding  at  the  contacts.  Ihe  tangential 
displacement  is  now  made  up  of  two  effects,  a  sliding  effect  and  a  con¬ 
tribution  due  to  slip.  On  the  other  hand,  the  norml  forces  and  dis¬ 
placements  are  still  connected  by  the  compliances. 

Furthermore,  since  Eq.  (5«44)  is  a  geometrical  relationship, 
it  holds  as  well  for  total  displacements  in  the  problem  at  band.  Hence, 

Ihe  vertical  strain  Eq.  (5.?^)  is  also  geometrical  and  ccm  be 
used  to  obtain  the  strain  for  the  medium. 


de 


zz 


C  dH 
2  2 


(3.46) 


Hext  we  note  that  the  equilibrium  Eqs.  (5.51a)  and  (5.31b), 
which  were  originally  written  in  terms  of  differential  stresses,  also 
hold  for  total  stresses.  Hence 
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/2 


zz 


(3.47a) 


'J  2 

«2  ■  ^2  -  T  ° 


(3.47t) 


We  have  already  determined  that  and  are  connected  by  the  coefficient 
of  friction  throughout  the  entire  loading,  i.e. 

(3.47c) 

Hie  behavior  of  a  granular  medium  of  face  centered  packing 
subjected  to  one  dimensional  compression  can  now  be  described  by  Eqs. 
(3.44),  (3.54),  (3.46),  and  (3-47). 

Equations  (3.47)  may  be  combined  to  eliminate  resulting  in 


(3.48a) 

Hjd-f)  -  '^2 

(3.48b) 

In  theoretical  and  applied  soil  mechanics  a  quantity  of  major 
Interest  in  one  dimensional  compression  is  K^,  the  ratio  between  the  hori¬ 
zontal  and.  vertical  stresses  at  rest.  From  Eqs.  (3*46)  this  ratio  is 
clearly 


K 


o 


1 

2 


(^ 

'l+f 


(3.49) 


Some  indication  of  the  variation  of  K  with  the  coefficient  of  friction 

o 

f  is  shown  in  Table  3.2. 

Let  us  now  determine  the  strains  resulting  from  one  dimensional 
compression.  Hie  substitution  of  Eq.  (5.2)  into  Eq.  (3.46)  yields 
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RN 

- 5m~ 


dN^ 


(5.50) 


Combining  Eqs.  (3-50)  and  (3*48a)  and  noting  from  Fig.  3.9  that 


0 


ZZ 


(3.51) 


gives 


^  r  2/3  _  ^-1/3 


(3.52) 


Integration  of  Eq.  (3-52)  yields 

_  2/3  2/3 

- 


ZZ 


zz 


(3.53) 


An  interesting  result  of  Eq.  (3*53)  is  that  the  stress-strain 
behavior  is  independent  of  the  radii  of  the  particles.  A  comparison  of 
the  stress -strain  curves  for  hydrostatic  and  one  dimensional  compression 
shows  that  the  curves  are  similar  in  shape,  but  turn  up  at  different 
rates.  Ihe  ratio  of  the  one  dimensional  strain  to  the  hydrostatic  strain 
at  the  same  level  of  stress,  a  ,  is  a  ratio  of  Eqs.  (3-53)  to  (3»20), 

ZZ 

\dilch  reduces  to 


(3-54) 


Valuett_of  f  of  0.1  and  0.2  give  values  for  of  1.88  and  1.77,  res¬ 

pectively.  lliis  indicates  that  the  hydrostatic  and  one  dimensional  stress- 
strain  curves  are  related  as  shown  in  Fig.  3*12. 
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5.2.2  Ifaloading  Cycle 

When  the  medium  is  loaded  the  tangential  forces  tend  to  resist 
the  sliding  motion.  For  unloading  the  tangential  forces  tend  to  resist 
the  sliding  movement  which  is  now  in  the  opposite  direction.  Hence,  the 
tangential  forces  reverse  their  direction  on  unloeuling.  Ihese  two  cases 
are  Illustrated  schematically  in  Fig.  5.13.  The  normal  and  tangential 
forces  are  still  related  by  the  coefficient  of  friction;  only  now  the 
direction  of  the  tangential  force  is  changed. 

The  new  equilibrium  equations  become 


N 


2 


N...  +  T., 


H 


(3.55) 


where 


(3.56) 


Consider  the  medium  to  be  loaded  from  0  to  fig*  3-1^ 

according  to  Eqs.  (5.53).  If  equilibrium  Eqs.  (3*55)  become  valid  at 

« 

the  instant  e  begins  to  decrease  from  €  ,  then  a  will  decrease  in- 

ZT.  zz  zz 

#  ' 

stantaneouBly  from  to  some  value  At  the  end  of  loading  and 

♦ 

0  are  related  by 
zz 


« 

0 

zz 


H 

2 


(3.57) 


\diereas  the  stress  o  is 
zz 


related  to  the  contact  forces  by  Eqs. 


(3.5*+) 


and  Eq.  (3.5l) 
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Ng  (1-^) 
'/2R^ 


(3.58) 


the  relation  between  a  and  a  at  the  Instant  \diere  e  Just  begins 

zz  zz  zz 

to  decrease  \rtien  the  value  has  essentially  the  same  value  while  the 
shearing  stress  has  switched  directions  Is 


t 

o 

zz 


*  /1-fx 

®zz  ^1+f^ 


(3.59) 


Combining  Eqs.  (3"55)^  (3*5l)>  (3"^)  *uid  (3-2)  yields  for  unloading 


(3.60) 


which  by  Integration  becomes 

2/3  2/3 

2  (3.61) 

Ihe  entire  stress -strain  history  for  one  cycle  of  loading  In 

,  « 

one  dimensional  compression  Is  shown  In  Fig.  3'^^  terms  of  o  and 

z  z 

1 

0  .  The  derivations  above  are  based  upon  the  assumption  that  the  tan- 

z  z 

gential  contact  forces  Immediately  reverse  directions  \rtien  unloading 
begins.  This  assumption  Is  not  quite  true,  however,  since  each  sphere 
must  exhibit  a  small  elastic  tangential  displacement  before  the  tangen¬ 
tial  forces  can  change  direction.  This  effect  however.  Is  small  with 
x"-3spect  to  the  tangential  displacement  due  to  sliding,  and  was  neglected 
in  the  analysis.  Due  to  the  above  assumption,  the  stress -strain  curve 
reilects  a  vertical  drop  In  stress  from  A  to  B  without  euiy  change  In 
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strain.  If  the  small  elastic  tangential  displacement  due  to  slip  were 
accounted  for,  the  stress -strain  curve  vould  not  have  a  cusp  in  it  at 
B,  and  would  follow  the  dotted  line  shown  in  Fig.  5.l4. 

Hie  stress -strain  behavior  shown  in  Fig.  3-14  exhibits  an 
energy  lose  after  a  cycle  of  loading  and  unloading  but  has  no  residual 
strain.  Hence,  this  medium  can  absorb  energy  without  any  permanent 
displacements. 

3.2.3  Energy  Absorption 

The  amount  of  energy  absorbed  by  the  medium  on  loading  and 
subsequent  unloading  can  now  be  determined. 

The  energy  required  to  load  up  to  an  applied  vertical  stress 
of  p^  and  the  strain  of  is 

"1 

^1  ^  yf  °  ZZ^^ZZ 
o 

(5.62) 

4  5/2 

Ej^  =  J  (1+f) 

Hie  energy  taken  out  of  the  medium  E^  during  unloading  back  to 
zero  from  a  stress  of  0^  and  strain  is 

Eg  =  I  (1-f)  (3.63) 


Hius,  the  energy  lost  is 


5/2 


AE  ■  E^  -  Eg  »  -  f€^ 


(3.64) 
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Likewise  the  ratio  of  the  energy  loss  to  the  energy  input  is 


^  _2f 

1+f 


(5.65) 


Uius,  a  very  significant  property  of  the  one  dimensional  stress -strain 
curve  is  that  the  per  cent  of  energy  absorbed  due  to  loading  and  unload¬ 
ing  is  always  constant  for  a  material  and  depends  only  on  the  coeffici¬ 
ent  of  friction  at  the  contacts.  Hence,  the  ratio  of  the  area  between 
the  loading  and  unloading  curves  to  the  area  under  the  loading  curve  is 
a  constant  given  by  Eq.  (3*65)*  per  cent  energy  absorbed  for  var¬ 

ious  coefficients  of  friction  is  given  in  Table  3.2. 

In  summary  it  should  be  clear  how  the  contact  forces  and  en¬ 


ergies  associated  with  this  medium  are  connected.  Bie  work  done  by  the 
normal  forces  during  deformation  s'tored  in  the  form  of  re¬ 

coverable  strain  energy.  On  the  other  hand,  the  work  done  by  the  tan¬ 
gential  forces  Tgdbg  is  a  nonrecoverable  energy  and  is  dissipated  as 
heat  into  the  medium.  As  seen  from  Table  3-2  this  energy  loss  during 


one  cycle  can  be  quite  significant.  In  fact,  with  a  coefficient  of 
friction  of  0.15,  the  dissipated  energy  is  26.1^  of  the  energy  put  into 


the  system. 

3.3  EquiveLLent  Discrete  Mass  Model  for  One  Dimensional  Static  and  Dynamic 


Behavior 

In  recent  years,  an  increasing  effort  has  been  devoted  to  study¬ 
ing  the  static  and  dynamic  behavior  of  soils  using  discrete  mass-spring 
models.  These  models  have  taken  many  shapes  with  various  contributions 
from  models  such  as  the  standard  Voigt  and  Maxwell  models.  These  model 
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studies  have  been  particiilarly  useful  in  wave  propagation  studies,  such 
as  the  work  of  Smith  and  Hewmark  I3.8].  Ihe  dynamic  equations  of  motion 
can  be  integrated  nvnnerically  with  a  digital  computer  using  the  S-method 
[3.9]  and  it  is  a  fairly  straightforward  approach  to  modify  soil  parame¬ 
ters  or  spring  stiffnesses  as  the  need  arises. 

As  might  be  expected  from  such  an  approach,  there  are  certain 
aspects  of  the  soil  behavior  which  escape  the  mass-spring  model.  Part 
of  this  inaccuracy  is  due  to  the  discreteness  of  the  system  and  part  can 
be  attributed  to  choice  of  the  model  itself. 

A  study  of  the  problems  associated  with  the  choice  of  various 
models  for  investigating  the  propagation  of  stress  waves  in  a  one  dimen¬ 
sional  medium  has  been  carried  out  by  Murtha  [3.IOJ. 

In  spite  of  the  Inherent  difficulties  associated  with  the  use 
of  models,  there  are  some  cases  where  the  model  is  the  only  hope  and  in 
fact  can  be  a  very  valuable  tool.  For  that  reason  the  authors  would 
like  to  suggest  a  one  dimensional  model  which  seems  to  exhibit  the  be¬ 
havior  expected  of  dry  coheslonless  sand,  and  is  based  upon  the  theory 
presented  in  Section  3-2. 

If  a  pressure  on  the  surface  of  the  earth  extends  over  a  large 
enough  area,  it  may  be  reasonable  to  consider  soil  completely  confined. 

If  this  is  valid  and  if  the  soil  composition  is  similar  to  sand,  the  one 
dimensioned  behavior  discussed  herein  may  very  well  be  a  good  approxima¬ 
tion  to  the  in  situ  soil  behavior. 

3.3>1  Horizontal  Model 

Figure  3- 15  gives  a  model  representation  of  the  sphere  system 
developed  in  Section  3-2.  Ihls  model  yields  the  same  stress -strain  cvur  e 
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ae  the  spheres  for  a  static  load  as  the  load  increases  or  decreases,  and 
the  equations  of  equilihrlum  in  the  limit  approach  those  of  the  spheres. 
The  model  consists  of  non-linear  spring  elements  and  Coulomb  damping  ele¬ 
ments  which  dissipate  energy  in  the  same  way  as  the  spheres.  The  reason 
that  it  depicts  horizontal  behavior  is  that  there  are  no  initial  stresses 
in  the  model  before  the  load  Is  applied.  Such  an  assumption  might  be 
reasonable  for  a  vertical  column  if  the  wei^t  of  the  soil  can  be  neglect¬ 
ed  compared  to  the  applied  stresses. 

3.3.2  Vertical  Model 

The  vertical  model  is  slightly  more  complex  them  the  horizontal 
one  in  that  the  weight  of  the  material  produces  initial  stresses  in  the 
model.  Since  the  stress -strain  curve  for  the  material  is  non-linear,  the 
stiffness  is  a  function  of  the  stress,  and  consequently,  a  function  of 
the  height  of  overburden.  The  overbvurden  pressure  increases  linearly- 
with  depth  and  it  can  therefore  be  Incorporated  in  the  equations  of  mo¬ 
tion. 

A  vertical  model  which  includes  the  initial  stresses  due  to 
the  overburden  of  the  material  is  given  in  Fig.  3 -16.  It  includes  a 
ciisnge  in  stiffness  associated  with  the  Increased  initial  stress  and  fol¬ 
lows  the  stress -strain  behavior  consistent  with  the  analysis  of  the 
sphere  medium  in  Section  3-2. 
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TABLE  3 . 1  VARIATION  OF  STRESS  AND  STRAIN  INVARIANTS 
(From  KJellman  [J.TJ) 


Av/v  =  (e^  +  Sg  +  e^)  = 

kg/  cm^ 

10-2 

10'^  in/ in 

1 

0.12 

0.04 

2 

0.24 

0.08 

4 

0.38 

0.13 

6 

0.47 

0.16 

6 

0.54 

0,18 

10 

0.62 

0.21 

12 

0.69 

0.23 

=  First  stress  Invariant 
=  First  strain  invariant 


TABLE  5.2  EFFECT  OF  THE  COEFFICIEMT  OF  FRICTION  ON 
ONE  DIMENSIONAL  STRESS-STRAIN  BEHAVIOR  OF 
FACE  CENTERED  ARRAY  OF  SPHERES 


f 

K 

0 

ae/e^  X  lOU- 

1  -  f 

2(l+f) 

g^xlOO 

0 

0.5 

0 

0.05 

0.45 

9.5 

0.10 

0.41 

18.2 

0.15 

0.37 

26.1 

0.20 

0.33 

33.3 

0.25 

0.30 

40.0 

0.50 

0.27 

46.2 

f  =  coefficient  of  friction  at  contact  points 
=  coefficient  of  earth  pressure  at  rest,  Pjj/p^ 

—  X  100^  =  precent  energy  absorbed  due  to  loading  and  unloading 
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Fig.  3.2  TWO  SPHERES  IN  NORMAL  CONTACT 


Fig.  3.3  DISTRIBUTIONS  OF  NORMAL  (d)  AND  TANGENTIAL  (t  ) 

COMPONENTS  OF  TRACTION  ON  THE  CONTACT  SURFACE 
OF  TWO  LIKE  SPHERES  SUBJECTED  TO  A  NORMAL  FORCE 
FOLLOWED  BY  A  MONOTONIC  TANGENTIAL  FORCE 
(after  Mindlin), 
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Fig.  3.  4  UNIT  CUBE  OF  A  FACE  CENTERED  CUBIC  ARRAY  OF 
EQUAL  SPHERES  SUBJECTED  TO  INCREMENT  FORCES 
AFTER  REF.  3.  3 


Fig.  3.  7  COMPONENTS  OF  FORCE  AT  A  CONTACT  AND  CORRESPONDING 

COMPONENTS  OF  DISPLACEMENT  OF  THE  CENTER  OF  A  SPHERE 
AFTER  REF.  3.  1 


kg. /cm. 
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Fig.  3.  8  VARIATION  OF  FIRST  STRESS  INVARIANT 
WITH  FIRST  STRAIN  INVARIANT 


z 


Fig.  3.  9  UNIT  CUBE  OF  A  FACE  CENTERED  CUBIC  ARRAY  OF  EQUAL 
SPHERES  SUBJECTED  TO  INCREMENTAL  FORCES  IN  ONE 
DIMENSIONAL  COMPRESSION 


Fig.  3.  11  VARIATION  OF  CONTACT  STRESSES  UNDER  ONE  DIMENSIONAL 
BEHAVIOR  WITH  AN  INITIAL  HYDROSTATIC  STRESS 


vertical  applied  stress, 
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Fig  3.  12  NON-DIMENSIONAL  STRESS  STRAIN  CURVES  FOR  A  FACE 

CENTERED  ARRAY  OF  SPHERES  SUBJECTED  TO  HYDROSTATIC 
AND  ONE  DIMENSIONAL  STATES  OF  STRESS 


(b)  Loading  Cycle  (c)  Unloading  Cycle 


Fig.  3.  13  SCHEMATIC  REPRESENTATION  OF  THE  DIRECTION  OF  THE 
CONTAC  T  FORCES  WHEN  LOADING  AND  UNLOADING  UNDER 
ONE  DIMENSIONAL  COMPRESSION 
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VERTICAL  STRAIN, 


Fig.  3.  14  LOADING  AND  UNLOADING  ONE  DIMENSIONAL  STRESS  STRAIN 
CURVE  FOR  A  FACE  CENTERED  ARRAY  OF  SPHERES 


.  15  ONE  DIMENSIONAL  HORIZONTAL  SAND  MODEL 
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chapter  4 

EFFECT  OF  GRAIN  CRUSHING  ON  THE  BEHAVIOR  OF  SANDS 

4.1  Introduction 

lyplcal  etress-strain  curvea  for  a  granular  material  such  as 
sand  subjected  to  static  loads  In  one  dimensional  compression  are  shovn 
in  Figs.  4.10  -  4.14.  The  general  shape  of  these  curves  has  been  quali¬ 
tatively  discussed  in  Chapter  1.  The  purpose  of  this  chapter  is  to  make 
some  evaluation  of  the  energy  absorbed  by  the  sand  in  the  pressure  ranges 
where  crushing  of  the  grains  takes  place.  One  should  note  that  this  is 
only  one  of  the  mechanisms  by  which  the  sand  absorbs  energy;  however,  en¬ 
ergy  absorption  in  this  region  can  be  very  significant. 

The  energy  absorbed  by  crushing  will  be  evaluated  by  determin¬ 
ing  the  new  surface  area  created  as  a  result  of  crxishing.  Ihe  increase 
in  new  surface  eurea  will  be  determined  from  a  statistical  evaluation  of 
the  change  in  the  grain  size  distribution  curve  for  each  pressure  range. 
Since  new  surface  area  can  be  related  to  energy  absorbed  in  the  crushing 
phenomena,  then  energy  absorption  due  to  the  creation  of  new  surfaces 
can  be  related  to  the  changes  in  the  grain  size  distribution  curve.  The 
results  of  this  analysis  led  to  additional  studies  of  stress-strain 
curves  for  sands  in  the  crushing  region,  >dilch  shoved  that  the  total  en¬ 
ergy  absorbed  due  to  particle  crushing  is  linearly  related  to  the  stress. 

4.2  Energy  Absorbed  Due  to  Crushing 

4.2.1  Statistical  Evaluation  of  Grain  Size  Distribution  Curve 

Since  the  change  in  the  grain  size  distribution  curve  for  a 
sand  is  a  measure  of  the  change  in  the  particle  sizes,  and  consequently, 
the  change  in  the  surface  area  of  the  particles,  it  is  also  a  measure 
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of  the  energy  absorbed  by  the  sand  due  to  crushing.  Since  each  group 
of  particles  Is  In  some  sense  a  random  accumulation  of  these  particles 
and,  furthermore,  since  the  crushing  of  the  particles  will  tend  to  oc¬ 
cur  in  a  somewhat  random  fashion.  It  would  appear  that  an  Investigation 
of  the  statlstlcsd  distribution  of  the  particles  and  certain  geometric 
parameters  Is  required. 

To  follow  this  thesis,  the  concepts  presented  by  Orr  and  Dalla- 
valle  [4.1,  4.2j  will  be  followed.  It  will  be  shown  how  certain  statis¬ 
tical  measurements  of  a  group  of  particles  may  be  used  to  obtain  useful 
parameters  ^Ich  characterize  the  distribution  of  these  particles. 

If  the  particles  are  all  spheres  (or  very  nearly  so)  the  aver¬ 
age  surface  area  S  may  be  computed  from  the  mean  linear  dimension  d  . 

s 

Hence 

® 

where 

f  z  frequency  with  idilch  a  particle  diameter  d  occurs 

n  B  total  number  of  particles  (n>>£f) 

Likewise,  the  mean  volume  V  is 


(4.2) 


If  due  allowance  is  made  for  weighting  of  size  axid  frequency, 
then  better  statistical  averages  would  be 


(4.3) 


I 
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V 


(4.4) 


\rtiere  d'  and  d'  are  denoted  the  mean  surface  and  mean  volume  diameters, 
respectively.  If  the  particles  are  not  spheres,  then  suitable  correc¬ 
tion  factors  called  shape  factors  must  be  used;  i.e.: 


S"  =  a  ( 5^)  =  a  d*'" 

8  \  n  /  SB 


(4.5) 


a  a  d-5 


v\  n  / 


V  V 


(4.6) 


Here  and  Ci^  are  the  surface -area  and  volume  shape  factors,  respec¬ 
tively.  Equations  (4.5)  and  (4.6)  are  the  correct  equations  to  use  if 
the  particles  are  irregularly  shaped. 

Some  typical  values  of  a  and  a  taken  from  [4.2j,  p.  27  are: 

6  V 


Cube  6  1 

Sphere  *  */6 

White  Seind  (smooth)  2. 1-2. 6 

Filter  Sand  (smooth)  2. 7-2. 9 

Crushed  Quartz  2. 1-2. 5  0.14-0.26 

In  general,  it  should  be  noted  that  while  there  is  a  somewhat  large  vcurla- 

tion  in  ct  and  cx  studies  have  shown  that  the  ratio  a  /a  is  usually  in 
s  v'  s  V  ^ 

the  range  of  6  to  7  and  closer  to  6. 

It  is  sometimes  of  interest  to  determine  the  specific  surface 
area  S  defined  as  the  surface  area  per  unit  volume  of  the  materlcd,  or 

V 


sometimes  the  surface  area  per  unit  wel^t  S^.  If  the  particles  are  all 
spheres,  these  turn  out  to  be 
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(4.7) 


i, 

Pd 


(4.8) 


vhere  P  is  the  density  and  d^  is  a  diameter  characteristic  of  specific 
surface  areas  as  shovn  in  Eq.  (4.l). 

By  investigating  Eqs.  (4.3)  and  (4.6)  it  is  seen  that  these 
equations  weigh  heavily  in  favor  of  larger  particles.  This  is  good  for 
volxime  or  mass  averages,  but  for  surface  areas  the  finest  particles  be¬ 
come  more  significant.  Hence  the  surface  area  computed  by  these  equa¬ 
tions  is  too  low.  To  overcome  these  effects  emd  to  obtain  a  better 
statistical  representation,  the  following  procedure  is  useful. 

Let  y  represent  the  volume  or  weight  fraction  of  all  particles 
measured  having  a  particular  diameter  d.  Then 


y  = 


(4.9) 


and 


(4.10) 


^ere  d^  Is  a  mean  weight  or  volume  diameter  of  the  distribution.  It 
should  be  noted  that  this  diameter  is  larger  than  that  calculated  from 
Eq.  (4.4),  since  the  contribution  of  the  larger  particles  Is  emphasized 
by  It. 

Nov  let  the  wel^^t  specific  surface  be  denoted  by  and  the 
specific  surface  be  denoted  by  S^.  Then 
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Lfd^S\ 

Efd^ 


From  Eq.  (4.8)  ve  have 


d  -  A.  _  6£fd^ 

p^w  pSfd^S’ 

V 

Since  =  6/pd  there  resvilts 


(4.11) 


(4.12) 


(^.13) 


The  usefulness  of  Eq.  (4.13)  is  that  It  eliminates  the  use  of 
shape  factors.  Equation  (4.13)  could  have  eulso  been  obtained  by  divid¬ 
ing  Eq.  (4.6)  by  Eq.  (4.5).  Hence  d  is  the  diameter  to  use  to  obtain 

8V 

the  surface  area  per  unit  volume  and  this  quantity  is  merely  6/d^^. 

Let  us  nov  consider  the  frequency  of  an  occurrence.  Ihxis  far 
we  have  been  dealing  with  measures  of  a  central  tendency  of  an  event 
happening;  however,  in  most  cases,  the  particle  size  distribution  will 
be  skewed  in  the  direction  of  increasing  size.  (Or  at  least  this  is 
what  is  to  be  expected  from  a  random  sampling  of  particles).  Ihese 
distributions  can  be  normalized,  however,  if  the  size  is  plotted  logari¬ 
thmically.  Thus  the  frequency  f  with  which  a  particle  of  diameter  d 
occurs  is: 


f 


Ln 


V2: 


In  0 


g 


exp 


In  d  -  In  M 


2  In^  o 

g 


(4.14) 


where 

M  ■  geometric  mean  diameter 
0  E  geometric  standard  deviation 

s 
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nxeae  are  defined  by 


in  M  =  E(i^) 


(4.15) 


In  a 


'(in  d  -  In  M)^ 
n 


(4.16) 


Afl  will  be  seen  later,  these  quantities  can  be  obtained  easily 
by  graphlcELl  means  from  the  grain  size  distribution  curve.  Furthermore, 
these  quantities  can  be  related  to  the  other  statistical  dimensions  as 
follows  where  the  logarithms  are  to  the  base  ten: 


log  d^^  «  log  +  4.605  log^  Og 
log  d'  ^  «  log  +  10.562  log^  0 

V  § 


log  dg^  -  log  M  +  5.757  log"  Og 


Also  if 


d  ■  mean  diameter  of  distribution 

m 

d^  ■  mean  weight  or  volvme  distribution 

there  results 


(4.17) 

(4.16) 

(4.19) 


d  -=  log  M  +  1.151  log^  0  (4.20) 

m  g 

d  -  log  M  +  8.059  log^  0  (4.21) 

The  above  equations  beglmlng  with  Eq.  (4.17)  are  called  the 


Hatch-Gboate  equations,  and  can  be  very  easily  applied  to  a  specific 
distribution  to  obtain  any  of  the  various  statistical  parameters  [4.5]. 
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4,2.2  Graphical  Determination  of  the  Geometric  Mean  Diameter  and 
Standard  Deviation 

If  the  diatrlbution  curve  for  a  group  of  particles  follows 
fairly  well  the  distribution  given  by  Eq.  (4.l4),  the  plot  on  a  log 
probability  grid  of  particle  size  versus  cumulative  percentages  (or 
fractions)  less  than  (or  greater  than)  a  specific  size  results  in  a 
stral^t  line. 

On  this  plot  the  size  corresponding  to  ^0^  Is  the  geometric 
mean  diameter  M.  Likewise  the  value  of  the  geometric  standard  devia¬ 
tion  a  can  be  related  to  the  sizes  at  84. 13^,  50^>  ojid  15* 8?)^  as  fol- 
fi 

lows : 


84.131^  size  50^  size 
®g  "  size  “  15.87^i  size 


(4.22) 


4.2.3  Particle  Size  Measxirements  Using  Sieves 

Hie  above  indicated  procedure  would  be  quite  simple  to  apply 
if  it  were  known  irtiat  the  number  of  particles  are  in  each  size  category. 
Gtafortunately,  this  is  not  the  case  with  sand  when  the  segregation  meth¬ 
od  is  carried  out  by  means  of  sieving.  When  sieving  is  used  the  results 
£ire  given  in  terms  of  percent  passing  (or  retained)  by  wel^t  and  not 
number . 

If  weight  percentages  are  used  cmd  the  plot  is  made  on  the 
log  probability  curve  the  result  will  still  plot  as  a  stralgiit  line  if 
the  distribution  follows  Eq.  (4.l4). 

On  this  weight  plot  we  define  the  median  and  standard  devia¬ 


tion  as 
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M'  ■  50^  size 

(4.23) 

,  84.13*  50*  size 

“  ~0r  “  i5.B7<  sire 

It  can  be  shown  [4.ll  that  these  values  are  related  to  M  and 


o'  ■  0 
g  g 

log  M  -  log  M'  -  6.908  log^  0  (4.24) 

g 

Hence  the  statistical  diameters  of  interest  can  be  obtained  as 


log  d^  -  log  M'  -  5.757  log^  Cg 

(4.25) 

log  d^  -  log  M'  -  3.454  log^  Og 

(4.26) 

log  d'  -  log  M'  -  4.605  log^  a 

*  g 

(4.27) 

log  d_^  -  log  M'  -  1.151  log^ 

g 

aunmary  it  should  be  noted  that  the 

(4.28) 

items  of  interest  in 

determining  only  the  surface  area  for  a  quantity  of  particles  do  not  in¬ 
clude  all  of  the  statistical  diameters  given  in  the  above  discussion. 
Since  is  a  meaawe  of  the  surface  area  per  unit  weight  it  is  the 
only  diameter  of  interest  for  this  study.  Hie  other  quantities  are  in¬ 
cluded  only  for  completeness. 

The  procedure  for  determining  the  surface  area  per  unit  volume 
of  a  material  is  to  make  a  weight  plot  on  a  log  probability  graph  and 
from  it  determine  the  median  and  standard  deviation  according  to  Eqs. 
(4.23).  From  this  information  d^^  can  be  obtained  using  Eq.  (4.26). 

Hie  surface  area  per  unit  volume  and  the  surface  area  per  unit  weight 


are  then 
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S 


V 


(4.29) 


(4.30) 


Since  Eq.  (4.29)  relates  to  the  volume  of  the  actual  solids, 
it  can  be  related  to  the  volume  of  the  loose  material  by  proper  use  of 
the  void  ratio  e  of  the  material. 

4.2.4  Energy  -  New  Surface  Area  Relationship 

In  recent  years  rather  extensive  studies  [4.4  -  4.12]  have  been 
carried  out  to  determine  the  energy  required  to  create  new  surface  area 
when  certain  solids  are  crushed.  Ihese  studies  have  been  carried  out  on 
such  materials  as  quartz,  glass,  florlte,  halite,  and  labradorlte.  The 
materials  were  subjected  to  both  static  loads  and  also  dynamic  effects 
created  by  dropping  weights  on  the  material.  Very  closely  controlled 
measurements  were  carried  out  to  determine  how  much  of  the  energy  input 
was  actually  used  in  the  crushing  of  the  solids. 

Considerations  were  given  to  temperature  changes,  deformations 
of  the  loading  apparatus  and  any  other  significant  energy  absorbing  mech¬ 
anisms.  In  general,  it  was  assumed  that  the  materials  did  not  undergo 
plastic  deformation  before  crushing,  but  only  exhibited  elastic  effects 
up  to  the  point  of  fracture.  It  was  felt  that  this  assumption  was  justi¬ 
fied  for  the  materlcds  noted  above,  since  under  normal  temperatiare  condi¬ 
tions  they  have  little  ductility. 


In  order  to  determine  the  surface  area  of  the  materials,  studies 
were  first  carried  out  using  methods  of  permeametry  as  discussed  in  l4.l]. 
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!Ihl6  Is  done  by  investigating  the  resistance  offered  to  a  flowing  fluid 
by  packed  particles  In  accordance  with  the  flow  lavs  of  D'Arcy.  In  some 
cases  the  fluid  used  vas  vater;  however,  for  smaller  particles  the  use 
of  air  in  place  of  a  fluid  Is  advantageous.  Ihe  permeability  method  for 
measuring  surface  area  has  the  advantage  of  being  a  relatively  simple 
procedure  to  carry  out;  however.  It  has  the  disadvantage  that  It  mea¬ 
sures  only  the  outside  surface  of  a  solid.  Such  things  as  cracks,  fis¬ 
sures,  and  pores  of  microscopic  and  submlcroscoplc  sizes  are  not  de¬ 
tected  by  this  method. 

In  view  of  the  above,  further  stisiles  in  this  series  used  gas 
adsorption  measxurementa  to  determine  the  surface  area  of  the  solids. 

This  technique  Is  also  discussed  In  [^.l]  and  Is  based  on  the  Idea  that 
the  surface  molecules  of  a  solid  are  bound  on  one  side  to  Inner  molecules 
but  are  Incompletely  attached  on  the  outside.  In  order  to  satisfy  the 
resulting  unbalance  of  atomic  and  molecular  forces,  the  surface  molecules 
attract  gas,  vapor,  or  liquid  molecules.  If  the  molecules  attracted  are 
those  of  a  gas,  the  phenomenon  Is  known  as  gas  adsorption.  It  Is  gen¬ 
erally  felt  that  gas  adsorption  measurements  provide  the  best  means  cur¬ 
rently  avallslble  for  determining  total  surface  area  of  a  solid. 

In  simple  terminology,  gas  adsorption  techniques  Involve  a  de¬ 
termination  of  the  quantity  of  a  gas  necessary  to  form  a  molecular  layer 
on  the  surface  to  be  measured.  The  number  of  molecules  required  to  form 
this  layer  may  be  evaluated,  azid  since  the  area  occupied  by  each  molecule 
Is  known  (or  may  be  estimated),  the  surface  area  of  the  material  may  be 
calculated. 

These  basic  techniques  have  been  applied  to  the  crushing  of 
single  particles  as  well  as  arrays  of  particles  for  some  of  the  materials 
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noted  above  [4.4  -  4,12].  Some  vai'latlons  were  made  In  the  loading  tech¬ 
niques  vlth  dynamic  as  veil  as  static  applications.  Tvo  graphs  summairlz- 
Ing  some  of  the  significant  data  related  to  the  behavior  of  quartz  were 
taken  from  [4.12]  and  are  given  as  Figs.  4.1  and  4,2. 

From  these  results  It  can  be  seen  that  the  new  area  per  work 
Input  was  measured  as  about  15.1  cmVkg  cm  In  single  particle  crushing 
compared  with  about  13.8  cm^/kg  cm  for  multiple  particle  crushing.  The 
higher  energy  required  to  crush  a  bed  of  particles  has  been  attributed 
to  the  friction  energy  losses  which  occur  when  a  group  of  particles  Is 
compressed. 

For  the  purpose  of  this  study  a  value  of  13.8  cm^kg  cm  or 
2. 46  In  /lb: In.  has  been  used  and  Is  considered  fairly  representative 
of  the  behavior  of  sand  as  well  as  quartz.  Unfortxmately,  there  are  no 
data  available  on  sand. 

4.2.5  Analysis  of  A  Grain  Size  Distribution  Curve  For  Sand 

The  concepts  discussed  In  the  previous  sections  have  been 
applied  to  determine  the  energy  absorbed  through  the  crushing  of  sand 
>rtien  subjected  to  large  pressures.  Unfortunately,  there  are  little  data 
available  In  the  llteratxure  which  can  be  utilized  for  correlation  pur¬ 
poses  and  to  date,  the  most  reliable  results  are  those  shown  In  Fig.  4.3 
which  were  taken  from  [4.13].  A  similar  set  of  data  Is  also  reported 
In  [4.14]. 

Fig.  4.3  gives  the  changes  In  the  grain  size  distribution  curves 
for  well-rounded  20-40  Ottawa  sand  when  subjected  to  various  ranges  of 
pressure  In  a  condition  of  one  dimensional  compression.  The  results  of 
an  analysis  of  these  data  are  given  In  Table  4.1  and  are  plotted  In  Fig. 
4.4. 
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The  results  In  Table  4.1  were  computed  by  the  following  pro¬ 
cedure.  In  general,  an  evalmtion  was  made  only  of  that  percent  of  mate¬ 
rial  'vrtiich  deviated  from  the  original  grain  size  distribution  curve. 

For  example,  at  the  1000  psi  range  only  about  3  percent  of  the  material 
was  affected,  ^ile  at  40,000  psi  about  88  percent  was  affected.  Even 
in  these  percentages  not  all  particles  were  crushed,  but  certain  statis¬ 
tical  average  diameters  were  changed  in  this  group.  An  investigation 
was  made  of  the  distribution  cxirve  of  these  altered  percentages  as  shown 
in  Figs.  4.5  -  4.8  and  the  specific  diameter  was  determined  for  this 
group  both  before  and  after  crushing  from  Eq.  (4.28),  From  these  diam¬ 
eters  a  determination  of  the  increase  in  surface  area  per  volume  was 
determined  from  Eq.  (4,29).  This  increase  was  converted  to  energy  from 
the  constant  2.46  in  /#  as  discussed  in  Section  4.2.4.  Ihese  calcula¬ 
tions  yield  the  energy  per  volume  of  solid  material  considered  in  the 
percentages  altered.  These  values  were  converted  back  to  energy  per 
total  volume  of  loose  material  assuming  the  void  ratio  to  be  0.60. 

It  is  shown  on  Fig.  4. 4  that  after  the  pressure  reaches  some 
2,000-3,000  psi,  the  energy  absorbed  due  to  crushing  seems  to  be  linearly 
related  to  the  pressure. 

Hiis  is  particularly  interesting  since  crushing  begins  to  have 
a  pronounced  effect  on  the  sand  behavior  at  about  this  pressure  even 
though  some  particles  begin  to  fail  in  the  region  of  5OO  to  1000  psi. 

From  these  meager  data  it  appears  that  the  rate  of  change  of  energy  to 
pressure  is  constant  and  is  approximately 


1  #in/ In^ 


0.01 


#Wln^ 

#/in^ 
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Uhfortiinately,  no  other  data  are  available  idilch  vUl  substantiate  this 
work.  Likewise^  the  grain  size  distribution  curves  in  Fig.  4.3  are  very 
likely  average  ctirves  for  many  tests.  Hence  this  study  must  be  viewed 
with  some  question  and  considered  as  indicating  only  a  qualitative  trend 
ixntll  further  experimental  results  can  be  obtained. 

4.3  Stress -Strain  Behavior  of  Sand  After  Crushing  Commences 
4.3.1  Experimental  Results 

The  statistical  study  carried  out  in  Section  4.2  Indicates 
that  in  the  pressure  region  after  grain  crushing  begins^  the  energy 
absorbed  by  sand  due  to  crushing  in  linearly  related  to  the  stress.  It 
was  pointed  out  that  this  energy  is  only  that  due  to  particle  crushing 
and  does  not  incluie  the  compaction  effect. 

Whether  these  two  pheno&ena  are  separable  effects  remains  to 
be  determined;  however,  they  are  distinctly  different  energy  absorbing 
mechanisms.  One  is  associated  with  the  breakage  of  the  moleculcur  bonds 
as  smaller  particles  are  formed  from  larger  ones.  The  other  effect  is 
the  energy  associated  with  the  settlement  and  compaction  due  to  a  change 
in  the  size  distribution  of  the  particles  which  results  in  an  Increased 
density. 

In  Section  4.2  it  was  shown  that  the  breakage  of  the  particles 
seems  to  follow  a  random  pattern.  If  it  is  assumed  that  the  newly  created 
fine  particles  are  transported  to  voids  in  the  medium  in  a  similar  ran¬ 
dom  pattern,  it  would  seem  plausible  that  the  total  energy  absorbed  due 
to  crushing  ml^t  also  be  linearly  related  to  the  stress. 

In  order  to  investigate  this  hypothesis,  let  us  study  some  ex¬ 
perimental  stress-strain  curves  for  sand  in  one  dimensional  compression. 
Figure  4.9  sh3vs  a  qualitative  stress-strain  curve  for  sand.  If  the  grains 
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of  sand  do  not  crush  the  stress-strain  curve  vUl  not  be  concave  dovn- 
vard  at  but  vlll  continue  to  be  concave  upward.  Since  crushing  does 
occur,  the  curve  changes  In  Its  shape.  If  the  loading  curve  Is  trans¬ 
lated  laterally  by  a  strain  of  e^  xmtll  It  cuts  the  crushed  stress -strain 
curve  at  a  stress  the  area  between  the  original  and  the  translated 
curve  Is  the  energy  absorbed  by  the  sand  due  to  crushing.  In  the  same 
way  the  energy  absorbed  due  to  crushing  at  some  stress  Is  the  cuaula- 
tlve  area  noted  In  Fig.  4.9. 

Figures  4.10  -  4.14  give  stress-strain  curves  from  some  tests 
on  20-40  Ottawa  Sand  carried  out  by  MIT  [4.13]  where  a  steel  chamber  was 
used  to  restrict  the  lateral  strains.  A  logarithmic  plot  of  these  same 
stress-strain  curves  In  Fig.  4.13  shows  that  prior  to  crushing,  the  stress 
is  related  to  some  power  of  the  strain  usually  between  1.3  and  1.6  and 
that  a  value  of  1.3  le  a  fairly  good  average.  The  theory  discussed  In 
Chapter  3  Tor  a  face  centered  array  of  spheres  predicts  this  power  to 
be  1.3. 

From  Fig.  4.3  It  is  an  easy  matter  to  extrapolate  the  stress - 
strain  behavior  ^Ich  would  occur  if  crushing  did  not  take  place.  Oils 
extrapolated  non-crushing  curve  Is  then  translated  along  the  strain  axis 
to  determine  where  It  crosses  the  true  curve  In  order  to  determine  the 
energy  absorbed  due  to  crushing. 

Figure  4.l6  shows  the  stress-strain  curve  of  Fig.  4.10  with  the 
translated  curves  extrapolated  from  the  pncrushlng  data.  The  energy  ab¬ 
sorbed  due  to  crushing  was  then  determined  by  planimeterlng  the  cumulative 
area  under  the  stress-strain  curve  corresponding  to  a  particular  stress. 

A  plot  of  the  energy  absorbed  due  to  crushing  Is  given  In 
Fig.  4.17  Tor  the  stress-strain  curve  In  Fig.  4.l6.  The  energy  absorbed 
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due  to  crushing  appears  to  be  linearly  related  to  the  pressure  for  a 
large  range  in  pressures.  Similar  studies  were  made  for  the  stress- 
strain  curves  given  as  Figs.  4.10  -  4.l4  and  a  plot  of  the  crushing  en¬ 
ergy  absorbed  versus  the  stress  is  also  given  in  Fig.  4.17.  There  is  a 
definite  trend  in  this  behavior  until  the  stresses  become  quite  large. 
Hence  it  can  be  concluded  from  the  meager  data  available  that  there  ap¬ 
pears  to  be  a  linear  relationship  between  stress  and  the  energy  absorbed 
due  to  crushing.  Further  study  needs  to  be  carried  out  on  this  work  as 
more  data  become  available;  however,  these  preliminary  investigations 
are  encovtraging. 

4.3,2  Theoretical  Stress-Strain  Curve  in  the  Crushing  Region 

As  was  discussed  in  the  previous  section,  there  is  experimental 
evidence  to  support  the  hypothesis  that  the  energy  absorbed  due  to  crush¬ 
ing  is  linearly  related  to  the  stress.  Slmllarlly,  a  study  of  the  stress- 
strain  curves  given  in  Figs.  4.10  -  4.14  as  well  as  the  experimental  re¬ 
sults  carried  out  on  this  study  and  presented  in  Chapter  3  indicate  that 
the  stress  prior  to  cirushlng  is  related  to  some  exponent  of  the  strain. 

Figure  4.18a  shows  a  qualitative  stress-strain  curve  for  sand 

with  a  discontinuity  at  o  the  crushing  stress.  This  stress  a  Is,  of 

c  c 

course,  not  a  well-defined  point  for  a  real  material;  however,  it  can  be 
determined  within  a  reasonable  range.  Let  the  stress  a  prior  to  crush¬ 
ing  be  expressed  as 

0  -  Oe^  (4.31) 

where  e  is  the  strain  and  a  and  3  are  experimented  constants. 

Likewise,  let  the  stress  after  crushing  begins  be  o^,  a  func¬ 
tion  of  the  strain.  Since  the  energy  absorbed  due  to  crushing  appears 
to  be  linearly  related  to  the  stress  we  have  from  Fig.  4.l8b. 
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dE 


7 


(4.32) 


vhere  E  la  the  energy  and  7  la  an  experlaental  eonatant.  But  the  Incre¬ 
ment  In  energy  dE  la  alao  the  Increment  In  area  under  the  atreaa-atraln 
curve  aa  ahovn  by  dA  on  Fig.  4.l6a.  Alao>  alnce  the  aldea  of  dA  are 
aaaumed  to  be  a  parallel 


dA  -  o^dc 


(4.33) 


Since  from  (4.32) 


7 


there  reaulta  from  (4.33)  And  (4.3>») 


(4.34) 


0^dC 


(4.35) 


or 


7d€  (4.36) 

®1 


Integrating  Bq.  (4.36)  ylelda  the  exponential  reault 


Aa 


1 


c 


7€ 


(4.37) 
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vhere  A  1b  a  constant  of  Integration.  Olils  constant  can  be  evaluated 

from  the  condition  that  at  the  crushing  stress  a  the  strain  Is  c  . 

c  c 

Hence  ve  obtain  flncdly 

=  o^e  “  ^c^  (*».38) 

Equation  (4.38)  satisfies  the  condition  that  the  energy  ab¬ 
sorbed  due  to  crushing  is  linearly  related  to  the  stress.  Tlie  constant 
7  can  be  obtained  from  curves  similar  to  those  given  In  Fig.  4.l8. 

Hence  analytical  expressions  for  a  one  dimensional  stress- 
strain  cxurve  take  the  following  form 

B 

0  =  Oe  prior  to  crushing 

(4.39) 

a  »  "  *c^  after  crushing 

Equation  (4.38)  Is  independent  of  the  stress -strain  behavior 
prior  to  crushing  and  Is  only  contingent  on  the  linear  energy  verstis  the 
Increment  of  stress  after  crushing.  Hence,  Eq.  (4.38)  may  not  be  lim¬ 
ited  to  only  one  dimensional  compression,  but  may  be  also  applicable  to 
other  states  of  stress. 

4. 4  Determination  of  the  Stresses  at  Which  Grain  Crushing  Occurs 

The  average  stress  for  crushing  depends  on  many  things  includ¬ 
ing  the  Initial  void  ratio  of  the  medium,  the  angularity  of  the  peurtl- 
cles,  the  duration  of  loading,  and  the  Inherent  strength  of  the  mineral 
which  composes  the  grains. 

For  the  one  dimensional  stress-strain  curve  of  a  given  material, 
the  effect  of  particle  breakage  appears  to  be  dependent  to  a  large  extent 
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on  the  initial  void  ratio  of  the  medium.  Ihls  statement  Is  substantiated 
in  Fig.  U.19  where  the  strains  of  Figs.  4.10  -  4.l4  have  been  "normalized" 
using  the  strain  at  20,000  psl  as  100^  strain.  These  data  show  that 
there  Is  a  general  trend  for  sand  with  the  higher  Initial  void  ratios  to 
crush  at  lower  levels  of  stress.  IRils  phenomenon  is  explainable  because 
at  a  high  initial  void  ratio  there  are  very  few  contacts  per  unit  volume, 
which  means  that  for  a  given  average  stress  the  contact  stresses  are 
hl£^er  in  a  sand  with  a  high  void  ratio  than  for  a  sand  with  a  low  void 
ratio.  Fui'ther  study  needs  to  be  csurried  out  on  the  determination  of  the 
crushing  stress  of  sand  and  hov  it  is  affected  by  particle  sizes  and  ar¬ 
rangements  . 
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Fig.  4.  1  WORK  INPUT  -  NEW  AREA  RELATIONSHIP  FOR  SINGLE  PARTICLE 
CRUSHING  OF  QUARTZ  BY  DROP- WEIGHT  AND  GLASS  BY  CALORI¬ 
METER  CRUSHER  (from  Zeleny  and  Piret  (4,  12)) 


Fig.  4.2  WORK  INPUT  -  NEW  AREA  RELATIONSHIP  FOR  MULTIPLE 

PARTICLE  SLOW  COMPRESSION  AND  DROP- WEIGHT  CRUSHING 
OF  QUARTZ  (fronr>  Zeleny  and  Piret  (4.  12) ) 


GRAIN  SIZE  DISTRIBUTION 
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GRAIN  SIZE  DISTRIBUTION  CURVES  FOR  OTTAWA  SAND 
(after  (4.  13)) 


Pressure,  1000  psi. 
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Energy  Absorbed,  Ib.-in/in? 


Fig.  4.  4  ENERGY  ABSORBED  DUE  TO  CRUSHING  OF 
20-40  OTTAWA  SAND 
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Fig.  4.5  LOG  PROBABILITY  PLOT  OF  GRAIN  SIZE 
DISTRIBUTION  CURVES 
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Fig.  4. 6 
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Fig.-  4.  8  log  probability  PLOT  OF  GRAIN  SIZE 
DIa  iRIBUTION  CURVES 


i'ig.  4.9  QUALITATIVE  STRESS  STRAIN  CURVE  FOR  SAND 
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Fig.  4.  10  '■TRESS  STRAIN  CURVE  FOR  SAND  IN 
ONE  DIMENSIONAL  COMPRESSION 
( from  data  in  (4.  13)) 


Stress,  cr,  1000  psi 
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Fig.  4.  U  STRESS  STRAIN  CURVE  FOR  SAND  IN 
ONE  DIMENSIONAL  COMPRESSION 
(from  data  in  (4.  13) ) 
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Fig.  4.  12  STRESS  STRAIN  CURVE  FOR  SAND  IN 
ONE  DIMENSIONAL  COMPRESSION 
(from  data  in  (4.  13)) 
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Fig.  4. 


13  STRESS  STRAIN  CURVE  FOR  SAND  IN 
ONE  DIMENSIONAL  COMPRESSION 
(from  data  in  (4.  13)) 


strain,  c,  Percent 


Fig  4.  14  STRESS  STRAIN  CURVE  FOR  SAND  IN 
ONE  DIMENSIONAL  COMPRESSION 
(from  data  in  (4.  13)) 
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Fig.  4.  16  cr- '  PARISON  OF  CRUSHING  AND  NON-CRUSHING 
BEHn-VIOR  OF  SAND 
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ENERGY  ABSORBED  DUE  TO  CRUSHING 
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.18  ASSUMED  STRESS  STRAIN  BEHAVIOR  FOR  SAND 


Percent  Of  Maximum  Strain,  Where 
«max.*  «  At  20,000  psi. 


Fig.  4.  19  NORMALIZED  STRESS  STRAIN  CURVES  FOR  OTTAWA  SAND 
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CHAPTER  5 

EXPERIMENTAL  PROGRAM 


5.1  Introduction 

An  experimental  program  for  the  purpose  of  measuring  the  pro¬ 
perties  of  real  soils  In  one  dimensional  compression  has  been  an  Integfral 
part  of  the  overall  efforts  of  this  research  study.  The  experimental  re¬ 
sults  from  these  tests  serve  several  purposes: 

1.  nie  basic  phenomena  observed  In  the  tests  can  be  compared 
with  the  conclusions  dravn  from  theoretical  analyses  of  Idealized  models 
^dilch  are  assumed  to  represent  soil. 

2.  The  energy  absorbing  capacity  of  soils  may  be  measured  by 
means  of  cyclic  loading. 

3<  Ibe  various  energy  absorbing  mechanisms  may  be  observed 
for  a  given  soil  at  different  pressure  levels  In  order  to  determine  the 
stress  ranges  vhere  each  mechanism  predominates. 

4.  Ibe  stress-strain  curves  obtained  may  be  used  to  study  one 
dimensional  vave  propagation  through  the  particular  soil  media  tested. 

3.  Ihe  lateral  stress  ^Ich  results  under  conditions  of 
zero  lateral  strain  can  be  measwed  vlth  the  «.>xperlmental  apparatus. 

6.  Ibe  effect  of  lateral  strain  on  the  principal  stress  ratio 

can  be  measured  by  allovlng  a  given  strain  In  the  horizontal  direction 
and  observing  the  change  In  vbere  Is  the  vertical  stress. 

7.  An  effective  Poisson's  ratio  for  soil  In  one  dimensional 
compression  can  be  evaluated  from  the  experimental  data  on 
relationship 


V 

l-» 


(5.1) 


A  device  was  designed  and  built  to  Investigate  the  one  dimen¬ 
sional  behavior  of  sand  in  the  high  pressure  regions.  The  apparatus 
can  measure  the  lateral  stress  under  conditions  of  zero  lateral  strain, 
and  provides  one  of  the  best  means  yet  developed  for  attaining  the 
condition  of  "zero"  lateral  strain.  Previous  Investigators  have  gener¬ 
ally  assumed  that  the  effects  of  small  lateral  strains  are  negligible, 
particularly  ^en  the  sample  Is  enclosed  In  steel  rings.  Research  by 
Speer  I5.1],  ^Ich  was  recently  pointed  out  by  two  of  the  authors  I5.2], 
shows  that  lateral  motion  significantly  affects  the  ratio  of  0g/oy 
The  results  of  Speers 's  research  is  presented  In  graphical  form  In  Fig. 
5.1  and  shows  that  a  lateral  displacement  of  4  x  10  ^  Inches  will  cause 
a  105t  reduction  in  the  value  of  Since  Speer's  work  was  with  a 

sand  sample  75/8  Inches  In  diameter,  the  above  diameter  change  corre¬ 
sponds  to  a  unit  strain  of  approximately  .5  x  lO"^  in. /in.  It  is  doubt¬ 
ful  to  the  writers  that  Speer  really  achieved  an  accxaracy  of  +  1  x  lO"^ 
in. /in.  since  small  temi)erature  fluctuations  would  cause  minor  varia¬ 
tions  of  at  least  that  magnitude  in  the  strain  gage  readings.  However, 
his  work  does  in  general  point  out  that  the  ratio  of  ®g/®Y  very  sensi¬ 
tive  to  lateral  movements.  This  phencxnena  had  eulready  been  observed  by 
Terzaghl  as  early  as  1934  In  connection  with  his  "Large  Retaining  Vail 
Tests"  at  MIT  [5.3].  Terza^l  concluded  that  an  outward  movement  of  the 
wall  of  .0007  h  In  the  case  of  a  well  compacted  dense  sand  was  enou£^ 
to  fully  mobilize  the  shear  strength  or.  In  other  words,  reduce  the 
ratio  of  ®g/ffy  “  ^  ^A  “  l+slny*  outward  movement  corresponds 

approximately  to  a  lateral  strain  of  about  12  x  lO”**^  in. /in.  Thus  for 
a  trxily  one  dimensional  test,  whereby  one  adso  wishes  to  measure  the 
magnitude  of  the  latersil  stresses  which  are  concomitant  with  the  vertical 
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load,  It  is  extremely  Important  to  closely  control  the  lateral  displace¬ 
ments,  In  fact,  preliminary  calculations  revealed  that  an  extremely 
thick  vailed  cylinder  subjected  to  an  internal  pressure  of  7,300  psl 
(vhich  is  an  approximation  to  the  lateral  soil  pressure  due  to  a  ver¬ 
tices  load  of  15^000  psi)  would  experience  a  radial  strain  of  5*25  x 
10  ^  in. /in.  if  the  specimen  were  7  inches  in  diameter  and  the  contain¬ 
ing  cylinder  was  assumed  to  have  an  infinite  external  radius.  Ihls 
strain  is  of  the  same  order  of  magnitude  as  the  tolerable  strains  listed 
by  Speer  and  Terzag^i.  llius,  in  order  to  study  one  dimensional  com¬ 
pression,  a  new  experimental  apparatus  had  to  be  designed  ^Ich  would 
restrict  the  lateral  deformations.  The  apparatus  developed  on  this 
program  is  discussed  in  the  following  sections. 

5.2  Experimental  Apparatus 

5.2.1  Description  of  the  Apparatus 

An  experimental  apparatus  va,B  designed  to  determine  the  stress- 
strain  relations  for  soil  under  one  dimensional  compression  and  to  mea¬ 
sure  the  lateral  stress  necessary  to  completely  restrain  the  sample.  Ihis 
apparatus  shown  schematically  in  Fig.  5*2  consists  essentially  of  a  thin 
steel  ring  idiich  contains  a  soil  sample.  Ihe  ring  is  surrounded  by  an 
annular  space  filled  with  oil  which  communicates  freely  with  hydraulic 
Jacks.  The  flexible  ring  and  oil  space  are  enclosed  in  a  thick  hollow 
cylinder  bolted  to  the  baseplate  in  order  to  withstand  the  high  fluid 
pressures. 

The  principle  upon  ^Ich  the  device  is  based  is  relatively 
simple.  As  the  vertical  load  is  applied  by  the  testing  head  there  are 
lateral  pressures  built  up  in  the  sand  which  tend  to  increase  the  diame¬ 
ter  of  the  thin  steel  ring.  Any  sllj^t  Increase  in  diameter  of  the  ring 
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Is  Immediately  Indicated  by  the  strain  gages  mounted  on  the  flexible 
ring  as  shovn  In  Figs.  ^.2  and  5*3 •  In  order  to  keep  the  lateral 
strains  zero,  the  oil  pressure  Is  modified  vlth  changes  In  the  verti¬ 
cal  load  In  such  a  manner  that  the  strain  Indicator  remains  balanced 
at  all  times  during  the  test.  When  the  strain  Indicator  remains  bal¬ 
anced  there  are  no  lateral  strains  and  the  oil  pressure  Is  equal  to 
the  lateral  soil  pressure  acting  against  the  side  of  the  container. 

The  apparatus  consists  of  the  following  fotir  basic  elements: 

(1)  A  thin  steel  ring  monitored  vlth  strain  gages  (Fig.  3*3a}> 

(2)  A  baseplate  (Part  B,  Fig.  5*2). 

(3)  A  thick  walled  cylinder  (Part  A,  Fig.  5*2). 

(4)  A  testing  head. 

The  testing  head  and  the  thin  steel  ring  contain  the  essential 
sensing  devices  for  making  the  desired  measurements  and  are  discussed  In 
the  following  sections. 

5.2.2  Fabrication  and  Caillbratlon  of  the  Strain  Gages  on  the  Steel 
Ring 

Four  Budd  Metalfllm  Strain  Gages  (Type  C6-II61)  were  mounted 
at  the  mid -height  of  the  steel  ring  at  Intervals  and  connected  In 
series  as  shovn  In  Fig.  5.3>  In  this  arrangement,  the  strains  sensed 
by  the  four  different  gages  are  averaged  since  the  change  In  resistance 
balanced  by  the  indicator  Is  the  sum  of  the  changes  In  all  four  gages. 

The  gages  are  "foil"  gages  with  a  gage  length  of  one  Inch  and 
a  grid  width  of  0.09  Inches.  Gages  of  this  proportion  were  chosen  In 
order  to  have  gages  relatively  Insensitive  to  axial  strains  In  the  ring 
arising  from  friction  between  the  soil  and  ring,  but  at  the  same  time 
very  sensitive  to  circumferential  strains  which  accompany  a  diameter 
change.  These  gages  worked  very  well  when  the  axis  of  the  one  Inch  gage 
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lengtb  vas  aligned  peirpendlcular  to  the  axis  of  the  ring,  lliese  gages 
were  also  very  convenient  because  they  are  extremely  thin  and  flat,  mak¬ 
ing  them  rather  Insensitive  to  the  all  around  pressure  of  the  oil  In 
which  they  are  submerged  during  test  conditions.  Nevertheless,  the 
effect  of  the  high  oil  pressure  on  the  gages  could  not  be  assumed  to 
be  negligible.  Ihe  gages  had  to  be  calibrated  for  these  effects. 

The  apparatus  Is  assembled  for  calibration  as  shown  In  Fig. 
^.4.  The  significant  featiire  of  this  assembly  Is  that  the  lower  "o" 
ring  between  the  bottom  of  the  steel  ring  and  the  baseplate  has  been 
removed  so  that  the  oil  In  the  annular  space  can  communicate  freely 
with  the  oil  Inside  of  the  sample  chamber.  A  steel  plug  one  Inch  thick 
Is  also  Inserted  Into  the  sample  chamber  to  confine  the  oil  and  a  test¬ 
ing  head  Is  lowered  flush  with  the  top  of  the  plug  to  supply  a  reaction 
of  sufficient  magnitude  to  keep  the  plug  In  place  during  calibration. 

Ihe  gages  were  calibrated  by  increasing  the  oil  pressure  In 
Increments  up  to  2,500  psl-  Since  the  "o"  ring  at  the  bottom  of  the 
steel  ring  was  omitted,  the  oil  pressure  In  the  annular  space  was  equal 
to  the  pressure  In  the  sample  chamber,  thus  giving  zero  net  pressure 
differential  across  the  ring  and  no  circumferential  stresses  or  strains. 
Hence  any  change  In  gage  reading  la  due  to  the  effects  of  the  all  around 
oil  pressure  on  the  gages.  Ihls  calibration  procedure  was  conducted 
several  times  and  the  calibration  curve  obtained  is  shown  in  Fig.  5*^* 
The  cujnre  was  reproducible  within  +  4  x  10~^  In. /in.  and  the  pressure 
effect  amounts  to  7  x  10~^  in./ln.  per  500  psl  of  oil  pressure.  Ccsnpen- 
satlon  for  these  pressure  effects  vas  made  when  the  tests  were  conducted 
on  sands. 

After  the  calibration  of  the  inner  diaphra^B  was  completed, 
the  device  vas  modified  for  testing  by  disassembling  the  apparatus  and 
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and  Inserting  the  "O"  ring  between  the  bottom  of  the  steel  ring  and  the 
baseplate  to  prevent  leakage  of  oil  i'rom  the  annular  space  into  the  sam¬ 
ple  chamber.  Hie  device  was  then  reassembled  as  shown  In  Fig.  ’p.2  for 
testing. 

5.2.3  Description  of  the  Testing  Head 

Hie  load  veis  applied  to  the  soil  sample  by  means  of  a  heavy, 
Internally  stiffened  piston,  mounted  In  a  120,000  lb.  Baldwin  hydraulic 
testing  machine  (Fig.  5.3b).  Hie  testing  head  Is  6.8OO  Inches  In  diame¬ 
ter  and  8.125  inches  high.  Hie  device  consists  of  two  rigid  steel  plates 
\dilch  are  welded  on  two  concentric  steel  cylinders  as  shown  In  Fig.  ^.6. 
Hils  figure  also  shows  the  manner  In  which  a  dynamometer  Is  Incorporated 
into  the  device  to  measure  the  pressure  over  the  center  square  inch  of 
the  loaded  area.  This  feature  enables  one  to  check  the  load  on  the  cen¬ 
ter  square  inch  against  the  average  load  over  the  entire  area.  Hie  pis¬ 
ton  and  the  dynamometer  were  designed  with  approximately  the  same  rela¬ 
tive  stiffness  in  order  to  maintain  a  uniform  deflection  of  the  specimen 
across  the  face  of  the  loading  device  as  the  sample  Is  compressed. 

The  vertical  displacement  of  the  sand  during  compression  was 
measxired  by  two  Ames  Dials  mounted  at  l80°  to  each  other  on  the  loading 
piston  as  shown  In  Fig.  5.7.  Hje  dials  measure  the  relative  displace¬ 
ment  between  the  movlhg  piston  and.  the  thick  vailed  cylinder  designated 
as  Part  A  in  Fig.  5.2.  Hiis  measured  relative  displiwement  Is  actually 
the  sum  of  the  vertical  displacement  of  the  soil  plus  the  strain  In  the 
testing  head  from  the  surface  of  loading  to  the  point  where  the  dials 
are  connected.  Hie  strains  in  the  head  were  so  small  compared  to  the 
strains  In  the  soil  that  they  were  neglected  and  the  measurement  was 
taken  to  represent  the  vertical  displacement  of  the  soil  sample.  Hie 


-109- 


dlals  aire  tu:eurate  to  l/ 10, 000  of  an  Inch  and  have  a  l"  travel.  VRien 
the  tvo  dials  read  differently,  the  average  reading  of  the  tvo  dials 
vas  taken  as  the  vertical  displacement  of  the  soil. 

5.2.4  Proportions  of  Test  Specimen 

The  proportions  of  the  test  specimen,  7  Inches  la  diameter  by 
2  Inches  hl^,  were  selected  to  minimize  the  effects  of  friction  between 
the  sand  and  the  steel  ring.  The  theoretical  basis  for  the  above  state¬ 
ment  Is  presented  In  the  following  analysis  from  l5-5]- 

Let  us  consider  a  confining  ring  of  radius  R  and  height  H  as 
shown  In  Fig.  5.8.  At  any  arbitrary  depth  z  below  the  surface  loaded 
by  the  force  P,  the  vertical  force  supported  by  the  soil  Is  designated 
eis  Q^.  The  force  may  be  expressed  In  terms  of  the  shearing  stress 
T_  and  the  applied  load  P  as: 

La 

z 


02  “  P 


2*Rt^z 


(5.2) 


The  shearing  stress  Is  sdso  related  to  by 
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f 
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\Aiere  A  Is  the  area  of  the  sample,  Is  the  ratio  of  the  horizontal  to 
vertical  stress  and  f  Is  the  coefficient  of  friction  between  the  soil 
and  the  ring. 

Substituting  Eq.  (5.3)  Into  Eq.  (5.2)  yields 
z 


f  dz 


(5.4) 
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Dlfferentlatlog  1x>th  sides  of  Eq.  (5*^)  and  rearranging  results  In 


-2*RK^f 


(5.5) 


Integrating  the  above  equation  between  the  appropriate  limits  ve  obtain; 


^  .  ,4 .  (,.5, 

p  o 

Rearranging  Eq.  (5.6)  we  obtain: 

-2tBK  fS  -2K  fH 

_ o  o 

A  R 

Qg  =  Pe  =  Pe  (5.7) 

-2K  fH 
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Equation  (5«7)  implies  that  If  e  R  approaches  1,  then  *  P,  and 
the  frictional  effects  become  negligible.  Uhls  relationship  shows  that 
the  r/h  ratio  should  be  as  large  as  possible  to  reduce  the  effect  of 
friction.  In  fact,  Increasing  the  r/h  ratio  by  a  factor  of  2  has  the 
same  effect  as  reducing  the  coefficient  of  friction  by  ^0^,  The  pre¬ 
sent  diameter  to  height  ratio  of  used  In  this  test  Is  considered 
sufficiently  large  to  minimize  the  frlctlonkl  effects,  but  on  future 
tests  It  would  be  desirable  to  vary  the  height  of  the  sample  in  order 
to  investigate  the  effect  of  this  ratio.  A  practical  limit  Is  i*esu:hed, 
however.  In  reducing  the  height  of  the  sample  since  the  hel£^t  of  the 
sample  Influencet  the  accuracy  of  the  vertical  strain  measurement.  If 
the  sample  becomes  too  thin,  the  Increments  In  the  displacements  can 
become  too  small  for  the  Ames  dials  to  sense,  especially  In  the  higher 
pressure  regions  ^ere  the  constrained  modulus  of  sand  approaches  ^00,000 
pel. 


5.?.1  General 

A  series  of  one  dimensional  compression  tests  with  lateral 
earth  pressure  measurement  was  conducted  on  a  dry  sand.  The  sand  was 
a  well  rounded)  coarse)  uniform)  silica  sand  from  Le  Suer)  Mlimesota. 

The  sand  is  very  similar  to  standard  Ottawa  sand  in  all  respects  except 
that  only  the  sizes  between  the  No.  10  and  No.  20  sieves  are  present 
in  this  material.  Relative  density  tests  on  this  material  showed  that 
the  void  ratios  for  the  loosest  and  densest  states  are  e  =»  .675 
and  “  *^55)  respectively.  The  loosest  state  was  obtained  by  poxir- 
ing  the  sand  in  a  container  of  known  volume  through  a  glass  funnel  to 
prevent  the  sand  from  falling.  The  densest  state  was  obtained  by  put¬ 
ting  the  sand  in  the  same  container  in  1/2"  layers  and  tapping  the 
sides  with  a  mallet  vmtil  it  appeared  as  if  the  surface  was  no  longer 
settling. 

The  initial  void  ratlO)  the  strain  rate  and  the  pressure  ranges 
of  cyclic  loading  were  varied  in  this  series  of  tests.  The  variation 
of  these  parameters  made  possible  a  stxidy  of  the  following  effects: 

(1)  The  influence  of  relative  density  on  K^. 

(2)  The  influence  of  strain  rate  on  K  . 

o 

(3)  The  relationship  between  the  overconsolidation  ratio  and 

(4)  The  effect  of  initial  void  ratio  on  the  stress -strain  curve. 

(5)  The  effect  of  initial  void  ratio  on  energy  dissipation. 

(6)  The  effect  of  strain  rate  on  the  stress-strain  curve. 

(7)  The  effect  of  pressure  level  on  energy  dissipation  due  to 
cyclic  loading. 


-112- 


5.3.2  Preparation  of  Teat  Speeinene 

Uben  preparing  the  test  specimens  it  is  extremely  important 

to  be  able  to  reproduce  the  same  initial  void  ratio  for  all  specimens, 

since  the  relative  density  significantly  affects  the  stress  deformation 

properties  of  sand.  Each  sample  of  sand  used  in  the  tests  consisted  of 

2,000  grams  of  oven-dry  material;  thus  the  volume  of  the  sand  solids 

vas  constant  for  each  test.  Since  the  void  ratio  e  is  by  definition 
V 

>diere  V  is  the  volume  of  the  voids  and  V  is  the  voltsae  of  the 
«  ^  V-V 

soil  solids,  it  may  also  be  written  —  %rtiere  V  represents  the  total 

s 

volume  of  the  sample.  Since  vas  constant  for  all  the  tests,  the 
void  ratio  vas  determined  by  controlling  the  total  volume,  V.  The  sam¬ 
ple  vas  cylindrlcaUy  shaped  so  this  vas  easily  accomplished  by  control¬ 
ling  the  height  of  the  sample. 

Each  test  specimen  vas  prepared  by  placing  2,000  grams  of 
oven-dry  sand  into  the  sample  chamber  and  Inserting  a  metal  plug  into 
the  chamber  above  the  sand.  This  steel  plug  fits  into  the  chamber  as 
shown  in  Fig.  5.^<  A  vibrator  vas  then  set  on  the  metal  plug  and  the 
sand  vas  vibrated  until  the  plug  settled  down  to  a  predetenslned  mark 
fixing  the  height  of  the  sample.  This  method  of  sample  preparation  has 
proved  to  be  rather  simple  and  convenient  and  the  test  results  thus  far 
indicate  that  very  good  reproducibility  has  been  achieved  by  the  employ¬ 
ment  of  this  technique. 

5.^,,  Test  Res\ilts 

5.4.1  Axial  Stress -Strain  Relationships 

The  axial  stress -strain  curves  for  each  of  the  eleven  tests 
are  shown  in  Figs.  5.9  through  5.19.  The  curves  show  the  entire  history 
of  loading  and  unloading  for  all  tests  Including  those  tests  where  the 
samples  were  subjected  to  seversd  ^eles  of  load. 
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nie  stress -strain  curves  In  generad  for  all  the  tests  are  non¬ 
linear  curves  which  are  concave  upward  on  both  the  loading  and  unloading 
portions.  The  non-llnearlty  of  the  curves  demonstrates  the  well  known 
fact  that  the  stiffness  of  a  given  sand  Is  highly  dependent  upon  stress 
level,  nils  stiffness  In  the  case  of  one  dimensional  compression  Is  usu¬ 
ally  designated  as  the  "Constrained  Modulus."  The  variation  of  the  con¬ 
strained  tangent  modulus  with  the  vertical  stress  for  the  Initial  loading 
In  tests  3;  ^  9>  end  13  are  shown  in  Figs.  ^.20  and  3*21.  These  curves 
also  show  the  effect  of  relative  density  on  the  constrained  modulus. 
Specimens  3  acd  U  had  an  Initial  void  ratio  of  0.62  and  had  a  constrained 
modulus  of  200,000  psl  at  a  stress  of  2,300  psi;  whereas  specimens  9  and 
13  had  an  Initial  void  ratio  of  0.^4  and  a  constrained  modulus  of  260,000 
at  2,300  pel  vertical  stress.  Figure  3 >22  shows  the  variation  of  con¬ 
strained  tangent  modulus  with  the  vertical  stress  for  Initial  loading  on 
all  tests  with  an  Initial  void  ratio  of  0.3^> 

The  stress-strain  curves  for  tests  3  and  4  are  Interesting  be¬ 
cause  In  a  stress  range  of  2,600  -  3,000  psl  the  initial  loading  curves 
show  a  tendency  for  the  curve  to  be  concave  downward.  Crushing  of  the 
grains  could  be  beard  at  this  stress  level  and  It  Is  possible  that  crush¬ 
ing  and  rearrangement  were  responsible  for  this  apparent  decrease  In  rig¬ 
idity.  This  tendency  toward  a  decreasing  modulus  In  the  2,800  psl  range 
was  not  observed  In  the  other  tests  even  though  the  crushing  could  be 
heard,  and  crushed  angular  particles  were  observed  after  testing.  The 
primary  reason  that  tests  3  ond  4  manifested  the  above  phenomenen  Is  be¬ 
cause  the  sample  for  these  tests  bad  an  Initial  void  ratio  of  0.62;  idiereas 
all  the  other  specimens  bad  an  Initial  void  ratio  of  0.34.  The  sand  with 
the  bigger  void  ratio  Is  free  to  rearrange  as  a  small  amovint  of  crushing 
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occurs,  and  the  strains  from  this  rearrangement  are  reflected  in  a  decreas¬ 
ing  modulus  as  in  tests  3  4.  The  sands  with  a  lower  void  ratio,  how¬ 

ever,  are  not  free  to  rearrange  to  a  great  degree  vAien  crushing  initiates 
and  thus  no  significant  strains  due  to  rearrangement  occur.  Even  the 
denser  sands  will  show  a  decreasing  modulus  at  some  stage  due  to  cjrushing, 
but  in  general,  one  would  expect  crushing  to  occur  at  a  higher  average 
stress  level  in  a  dense  sand  than  in  a  loose  sand.  A  dense  sand  has  more 
contact  points  per  unit  volume  than  a  loose  sand  and  thus  for  the  same 
average  applied  stress,  the  loose  sand  has  higher  grain  to  grain  contact 
stresses  than  the  dense  sand.  Therefore,  the  loose  sand  begins  crushing 
at  a  lower  stress  than  the  dense  sand.  The  phenomena  discxissed  above  is 
Illustrated  very  well  In  Fig.  4.19  of  this  report. 

The  loading  and  subsequent  unloading  of  the  sample  showed  that 
on  the  average  about  75/^  of  the  axial  strain  was  recoverable.  This  strain 
was  probably  the  result  of  elastic  deformations  at  the  points  of  contact. 
The  non -recoverable  portion  of  the  strains  most  definitely  was  due  to  ir¬ 
reversible  rearrangement  of  grains.  The  significant  portion  of  the  re¬ 
arrangement,  however,  takes  place  on  the  first  cycle  of  loeuiing  since  sub¬ 
sequent  cycles  are  nearly  reproducible  and  show  only  a  slight  tendency  to 
rearrange  further.  These  phenomena  are  illustrated  fairly  well  by  test  8. 

The  stress -strain  data  from  the  tests  are  presented  on  graphs  of 
log^Q  stress  versus  log^Q  strain  in  Figs.  5 -25  through  The  stress - 

strain  data  plot  is  a  straight  line  on  these  graphs  which  means  that  the 
axial  stress  may  be  expressed  in  terms  of  the  axial  strain  by 


(5.8) 
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The  above  expression  may  also  be  written  as: 

^°®10  “  ^°®10®Z 


(5.9) 


Thus  from  Eq.  (5.9)  it  is  apparent  that  the  parameter  n  is  the  recipro¬ 
cal  of  the  slope  as  measured  from  the  graphs  shown  and  K  is  the  value  of 
the  strain  at  which  the  stral^t  line  Intersects  the  strain  axis. 

The  values  of  n  range  from  O.58  to  0.53  for  these  tests  as  compared  to  the 
value  of  2/3  which  was  predicted  for  n  from  the  theoretical  analysis  pre¬ 
sented  in  Chapter  3-  The  results  do  show,  however,  that  an  exponential 
relationship  does  occur  as  predicted  by  theory,  but  the  value  of  n  from 
tests  seems  to  be  consistently  less  than  2/3. 

The  theory  from  Chapter  3  also  predicted  a  relationship  between 
the  constrained  tangent  modulus,  M  ,  and  the  vertical  stress,  a  ,  as 
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where  C^  is  a  constant.  Figure  5.3*+  shows  the  variation  of  log  and 
log  0^^  for  all  tests  on  samples  with  an  InitleJL  void  ratio  of  O.54. 
Since  the  relationship  is  a  straight  line  the  actual  relationship  be¬ 
tween  the  constrained  tangent  modulus  and  the  vertical  stress  is 

M  -  Cc  (5.11) 

c  zz 


where  k  is  the  slope  of  the  line.  The  value  of  k  for  this  set  of  experi¬ 
ments  is  1/2  whereas  the  theory  from  Chapter  5  predicted  that  k  should  be 

1/3. 

Ihe  small  number  of  tests  conducted  thus  far  are  not  sufficient 
to  make  any  definite  conclusions  about  the  effect  of  strain  rate  on  the 
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stresB -strain  curves;  however,  there  appears  to  be  a  definite  trend  even 
though  the  range  of  strain  rate  has  been  relatively  small.  Ibe  constrain¬ 
ed  moduli,  in  the  2,^00  -  3>000  psl  vertical  stress  range,  for  samples 
with  an  Initial  void  ratio  of  O.34  were  compared  for  loading  rates  vary¬ 
ing  from  0.003  to  0.0^  In. /min.  llie  results  are  shown  below. 


Ibe  results  of  the  lateral  earth  pressure  measurements  for  the 
Initial  loading  of  each  specimen  are  shown  In  Figs.  3*3^  through  ^.43. 
The  lateral  pressure  Is  plotted  as  ordinate  and  the  vertical  stress  as 
abscissa.  The  experimented,  data  plot  In  a  strai^t  line  for  all  curves 
and  K^,  the  coefficient  of  earth  pressure  at  rest,  is  numerically  equal 
to  the  slope  of  the  stralc^t  line. 

The  Influence  of  relative  density  on  can  be  seen  by  compar¬ 
ing  for  the  samples  having  an  inltled  void  ratio  of  0.62  with  the 
sample  having  an  inltled  void  ratio  of  O.34.  The  samples  with  the  hlc^ 
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void  ratio  had  measured  values  of  ranging  from  0.390  to  0.415;  vhereaa 
the  denser  samples  tested  at  the  same  rate  of  loading  as  the  above  samples 
had  values  from  0.540  to  0,360.  Thus  appears  to  decrease  as  the 
relative  density  of  the  sand  Increases,  nils  veorlatlon  vas  expected  since 
It  Is  veil  knovn  that  the  angle  of  internal  friction  Increases  as  the 
relative  density  Increases.  Other  Investigators  such  as  Jaky  [5.4]  have 
shown  that  the  value  of  decreases  as  the  angle  of  Internal  friction 
Increases  for  a  soil  vlth  an  over -consolidation  ratio  of  1. 

The  value  of  K  for  a  soil  of  constant  Initial  void  ratio  also 
o 

varies  slightly  vlth  strain  rate.  The  general  trend  of  the  data  shove 
that  decreases  vlth  Increasing  strain  rate.  Tests  on  the  dense  sand 
(e  »  0.54)  shoved  that  varied  from  0.375  for  a  consolidation  type 
test  as  In  test  10  to  as  low  as  0.350  for  the  fastest  loading  rate  of 
0.004  In. /min.  The  values  of  for  the  Intermediate  loading  rates  of 
0.01  and  0.02  in. /min.  were  0.350  and  0.355>  respectively. 

Thus  far;  the  discussion  of  the  coefficient  of  lateral  earth 
pressure  at  rest  has  been  restricted  to  that  portion  of  the  stress -strain 
curve  where  the  ratio  of  the  maximum  previous  stress  to  the  existing  stress 
on  the  sand  Is  unity  (l.e.,  previous  stress  _  Oygrcongoiidution  ratio  * 
OCR  «  1).  In  this  region  the  coefficient  of  earth  pressure  at  rest  Is  a 
constant;  as  predicted  by  the  theory  of  Chapter  5;  but  If  the  specimen  is 
unloaded  from  some  maximum  pressure;  then  the  OCR  becomes  greater  than  1 
and  the  value  of  does  not  remain  constant  as  unloading  progresses.  The 
variation  of  with  XR  Is  illustrated  In  Fig.  5*46  for  tests  3  and  4. 

This  graph  clearly  shovs  that  Increases  vlth  Increasing  OCR.  The  maxi¬ 
mum  value  of  measured  during  unloading  vas  2.0  and  It  appeared  to  be 
still  Increasing.  If  measurements  could  be  more  accurately  obtained  In 


the  final  stages  of  unloading  ^ere  the  pressures  are  very  lov,  It  would 
probably  be  found  that  the  values  of  would  ultimately  approach  the 
coefficient  of  passive  earth  pressure. 

Figures  5*^7  and  ^.kQ  show  the  variation  of  lateral  earth  pres¬ 
sure  with  vertical  stress  for  1  ij 2  cycles  of  loading  In  tests  12  and  13 . 
These  curves  show  that  with  cyclic  loading  the  coefficient  of  lateral 
earth  pressure  can  gradually  be  built  up  with  each  cycle  of  load.  Ob¬ 
serve  the  data  for  test  13  for  example  at  points  a  and  b.  The  value  of-*- 
at  a  Is  higher  than  at  b  even  though  the  overconsolldatlon  ratio  Is 
the  same  at  these  two  points.  Thus,  the  value  of  for  a  given  sand  is 
not  only  a  function  of  stress  history  as  expressed  by  the  OCR,  but  it 
also  depends  on  the  number  of  preceding  cycles.  Further  tests  are  ex¬ 
pected  to  indicate  how  much  can  be  Increased  at  a  given  OCR  by  cyclic 
loading. 

5.4,3  Energy  Absorption 

The  energy  absorption  characteristic  of  the  dry  sand  were  In¬ 
vestigated  by  means  of  cyclic  loading.  In  general  the  first  cycle  was 
conducted  by  loading  to  a  maximum  pressure  of  3^290  psl  and  iinloadlng 
to  zero.  The  deformations  in  the  first  cycle  were  only  about  751^  re¬ 
coverable  due  to  the  fact  that  rearrangement  of  the  grains  caused  ir¬ 
recoverable  strains.  The  sample  was  then  loaded  to  some  lower  stress 
level,  usually  1,100  psl,  and  then  unloaded.  It  was  found  generally  that 
the  second  cycle  deformations  were  almost  100^  recoverable,  but  energy 
was  still  being  lost  because  the  loading  and  unloading  paths  were  differ¬ 
ent.  If  the  second  cycle  was  followed  by  a  third  cycle  Identical  to  the 
second  (such  as  In  test  3,  5^  6,  "J,  and  6),  It  wu  found  that  the  third 
cycle  traversed  essentially  the  sane  loop  with  very  little  Irrecoverable 
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straic  but  that  a  significant  antount  of  energy  vets  still  being  dissipated 
due  to  the  difference  between  the  loading  and  unloading  paths.  In  fact^ 
these  loops  were  so  nearly  Identical  that  for  tests  ^  and  7  the  sec¬ 
ond  and  third  cycles  traversed  exactly  the  same  path  and  the  small  loop 
shown  on  the  stress -strain  curves  for  these  tests  actually  represents 
both  the  second  and  third  cycles. 

Ihe  energy  lost  in  hysteresis  loops  for  various  tests  was 
evaluated  by  computing  the  area  enclosed  In  the  loops  traversed  In  the 
second  and  third  cycles  of  loading.  It  was  found  convenient  to  express 
the  energy  lost  as  a  percentage  of  the  energy  input.  Hie  results  of 
these  measurements  are  shown  In  Table  5.1*  iQ  general  It  was  found  that 
the  ratio  of  energy  lost  to  energy  Input  Is  practically  a  constant  re¬ 
gardless  of  the  magnitude  of  vertical  stress,  and  varies  between  the  nar¬ 
row  limits  of  0.275  ^  0.33^.  For  samples  irtilch  had  the  same  history 
of  loading,  It  was  found  that  the  samples  with  the  higher  Initial  void 
ratios  always  lost  more  energy  than  those  with  the  lower  initial  void 
ratio . 

The  fact  that  the  energy  lost/  energy  Input  ratio  has  been  ob¬ 
served  to  be  fairly  constant  regardless  of  the  range  of  stress  is  signi¬ 
ficant,  but  this  relationship  Is  restricted  to  the  maximum  value  of  ver¬ 
tical  stress  employed  In  these  tests.  At  higher  stress  levels  crushing 
may  be  a  major  factor  and  thus  other  mechanisms  of  energy  dissipation 
enter  Into  the  picture.  Further  tests  Into  the  higher  pressure  regions 
will  enable  a  better  determination  of  the  range  of  pressures  for  idilch 
this  energy  ratio  Is  applicable. 

The  theory  of  Chapter  3  predicted  that  the  ratio  of  energy  ab¬ 
sorbed  to  Input  energy  for  a  hysteresis  loop  from  a  stress  of  zero  to  any 
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stress  belov  the  crushing  stress  should  be  a  constant  related  to  the  co¬ 
efficient  of  friction  between  grains.  Ibus  the  actual  meuurements  are 
in  agreement  with  theory  because  this  ratio  was  observed  to  be  practically 
a  constant.  The  laboratory  measurements,  however,  provide  another  check 
on  the  theory  iirtilch  is  quite  interesting.  Tie  theory  predicted  that 


and 


■'o '  <Er> 


(5.12) 


^  _2f 

E  “  1+f 
1 


(5.13) 


Since  K  and  ^  have  both  been  measured  in  the  laboratory  tests, 

O  * 

Che  coefficient  of  friction  or  more  appropriately,  the  pseudo  coefficient 
of  friction  can  be  evaluated  by  two  independent  calculations  for  each  test. 
If  the  theory  is  reasonably  appropriate  the  two  values  of  f  should  be  in 
good  agreement  for  each  test.  Ihese  calculations  were  made  for  each  test 
and  are  presented  in  Table  5.2.  Tests  5  and  4  do  not  show  good  agreement 
for  the  two  calculated  values  of  the  coefficient  of  friction.  Tests  5,  6, 
7,  8,  9,  11,  12,  and  13,  however,  show  a  very  good  correlation  between  the 
two  Independently  calculated  values  of  the  coefficient  of  friction.  The 
latter  group  of  tests  had  an  Initial  void  ratio  of  0.54  as  compared  to  0.62 
for  tests  3  and  4.  Since  the  theory  of  Chapter  3  assumes  a  granular  medium 
in  a  dense  state,  it  is  not  surprising  that  the  tests  on  sand  in  the  denser 
state  shovdd  give  better  correlation  with  the  theory  than  tests  as  higher 
initial  void  ratios  idiich  allows  more  rearrangement  of  the  grains. 
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TABLE  5.1 

ENERGY  LOSSES  FOR  REPEATED  LOADINGS 


Test  No. 

Loop  Range 
psi 

Mean  Stress 
psi 

Energy  Lost 

> 

3 

0  -  1100 

548 

33.4 

4 

0  -  1100 

548 

33.1 

5 

0  -  1100 

548 

32.8 

6 

0  -  1100 

548 

30.2 

7 

0  -  1100 

548 

27.9 

8 

0  -  1100 

548 

27.5 

8 

0  -  32190 

i64t 

21.3 

9 

0  -  1100 

vn 

CO 

28.7 

11 

0  -  2195 

1099 

28.6 
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TABLE  5.2 

CORRELATION  OF  t,  K^,  AND  AE/e^ 

1-2K 
t  = 


3 

.62 

.12 

.20 

4 

.62 

.09 

.20 

5 

.54 

.18 

.20 

6 

.54 

.16 

.18 

7 

.20 

.16 

8 

.20 

.16 

9 

.54 

.19 

.17 

10 

.54  (No  Hysteresis  loop  on  this  test) 

11 

.54 

.18 

.17 

12 

.54 

.18 

.14 

13 

.54 

.18 

.15 

1/2  E^ 

(1  -  1/2  f) 
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Fig.  5.5  CALIBRATION  CURVE  FOR  THE  EFFECT  OF  OIL  PRESSURE 
ON  THE  METALFILM  STRAIN  GASES 
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CROSS-SECTION  OF  TESTING  HEAD 
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Fig.  5.8  SCHEMATIC  DIAGRAM  ILLUSTRATING  THE  EFFECT  OF 
SIDE  FRICTION  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  psi. 
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Fig.  5.9  STRESS  STRAIN  CURVE  FOR  MINNESOTA  SAND 
IN  ONE  DIMENSIONAL  COMPRESSION 
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Fig.  5.  n  STRESS  STRAIN  CURVE  FOR  MINNESOTA  SAND 
IN  ONE  DIMENSIONAL  COMPRESSION 
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Fig.  5.  17  STRESS  STRAIN  CURVE  FOR  MINNESOTA  SAND 
IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  o-w,  psi. 
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Fig.  5.  21  CONSTRAINED  MODULUS  VS.  VERTICAL  STRESS 
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,  23  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  psi. 
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Fig..  5.  24  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINMESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  a^,  psi. 


Fig.  5.  25  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  psi. 


•149- 


Fig.  5.  26  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  o-^,  psi. 
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Fig.  5.27  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  cr^,  psi. 
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Fig.  5.  28  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  ffw,  psi. 
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Fig.  5.  29  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  o-w,  psi. 
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Fig.  5.  30  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  o-w,  psi 
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Fig,  5.  31  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Stress,  o-w,  psi. 
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Fig.  5.  32  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  Strain,  Cy,  in. /in. 


Fig.  5.  33  LOGARITHMIC  PLOT  OF  STRESS  STRAIN  CURVE  FOR 

MINNESOTA  SAND  IN  ONE  DIMENSIONAL  COMPRESSION 


Vertical  ttress,  pal 


Fig*  5*34  COHSTRAUIED  MODOLUS  VS  VERIK 
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Fig.  5.40  HORIZONTAL  STRESS  VS.  VERTICAL  STRESS  FOR  INITIAL 
LOADING  OF  MINNESOTA  SAND  IN  ONE  DIMENSIONAL 
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PRESSURE  AT  REST,  K  ,  AND  THE  OVERCONSOLIDATION 
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Fig.  5.47  HORIZONTAL  STRESS  VS.  VERTICAL  STRESS  FOR  ONE 
AND  ONE -HALF  CYCLES  OF  LOADING 
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APPENDIX  A 

APPLICATION  OF  THE  HERTZ  THEOKf  TO 
THE  BEHAVIOR  OF  A  GRANULAR  MEDIUM 

A  question  of  paramount  Importance  aurlsea  In  the  development  of  a  . 
theory  of  granular  media  based  on  the  Hertz  theory  [a.i]  .  It  must  first  be 
concluded  that  the  Hertz  theory  Is  valid  In  the  pressure  regions  of  Interest 
before  the  theory  has  any  usefuUness.  In  vlev  of  this  it  may  be  vorthvhlle 
to  review  the  developnent  of  the  Hertz  theory  here  for  the  convenience  of  the 
reader  aiod  note  slmultaneoiuly  the  fundamental  assimptlorrs  as  they  arise. 

The  derivation  given  in  the  following  is  a  siaaiary  of  the  presentation 
given  by  Timoshenko  and  Goodler  [a.2J^  and  the  asswptlons  are  discussed  by 
the  writers  In  relation  to  applying  the  results  to  a  granular  medlua  composed 
of  well  rounded  quartz  sand. 

Let  us  first  look  at  the  pressure  between  two  spherical  bodies  held 
In  contact  by  a  nozmal  force.  In  the  solution  of  this  problem  it  is  assumed  that 
at  the  point  of  contact  these  bodies  have  spherical  surfaces  with  the  radii  R^ 
and  Rg  (Fig.  A.l).  If  there  is  no  pressure  between  the  bodies  we  have  contact 
at  one  point  0.  The  normal  distcuxces  from  the  tangent  plane  at  0  to  points  ' 
such  M  M  and  N,  on  a  meridian  section  of  the  spheres  at  a  very  small  distance 
r  from  the  axis  z^  and  can  be  approximated  in  the  following  wanner  (Fig.  A.2): 

z  «  u  tan  ^/2  ~  1/2  u  tan  B 

z  -  1/2  u  .  u/B  ■  u^/2R  (A.1) 

z  ■  r^/2R  where  r®  -  u® 


Thus  in  (Fig.  A.l) 


173 


*1  “  2Rj^  ^  *2  “  aig 

ISm  dlstaaee  b«tiM«n  points  N  and  M  is 

r2(R^  +  R^) 


(A.2) 


(A.3) 


In  the  particular  ease  of  oontaet  betveen  two  spheres  of  equal  radius, 
R,  we  hare: 


*1+  *2 


2 

r^  2H 


2R‘ 


(AA) 


If  the  bodies  are  pressed  together  along  the  nomal  at  0  by  a  force 
P,  there  will  be  a  loeal  defozmtlon  near  the  point  of  contact  producing  contact 
orer  a  saall  surface  with  a  circular  boundary,  called  the  surface  of  contact. 
Assmlng  that  the  radii  of  eurrature  R^  and  R^  are  rery  large  in  conparison  with 
the  radius  of  the  boundary  of  the  surface  of  contact, we  can  apply,  in  discussing 
local  deforaation,  the  results  obtsdned  for  a  point  load  on  a  seai->lnflnite 
boundary.  Let  denote  the  displacasent  due  to  the  loeal  defoxaation  in  the 
direction  of  a  point  such  as  N  on  the  surface  of  the  lower  ball  (Pig.  A.l), 
and  Vg  denote  the  sesw  dXsplaceaent  in  the  direction  for  a  point  stich  as  N 
of  the  vppwr  ball.  If  it  is  asstseed  that  the  tangent  plane  at  0  reaalns  laaovable 
dxtring  loeal  compression,  then,  due  to  this  eoapression,  asQr  two  points  of  the 
bodies  on  the  axes  s^  and  at  large  distances*  ftrosi  0  will  approach  each  other 


*Sueh  distances  that  defomations  doe  to  the  eoBvression  at  these  points  can 
be  negleeted. 


Vjr  a  eartaln  cnnount  a,  and  the  dletaace  between  two  points  such  as  N  and  H 
will  diffllnlsb  by  a  -  (v^  -f-  v^).  If  the  two  spheres  are  of  equal  radius, 
the  distance  between  such  points  as  N  and  H  will  dlsilnish  by  a  -  2w.  If 
finally,  due  to  local  coapression,  the  points  M  and  H  coaw  inside  the  surface 


of  contact,  we  have 


+  Vg) 


2l+  Zg 


+  1^2^ 


2«1«2 


(A.5) 


If  the  scares  have  a  radius  R  >  then  we  have 

a  -  2w  -  2z  »  r^R 


Ihus 


(A.6) 


at  any  point  r  on  the  surface  of  contact. 

Let  us  now  consider  local  deformations.  From  the  condition  of  syanetry 


it  can  be  concluded  that  the  Intensity  of  pressure  q  between  the  bodies  in  contact 
and  the  corresponding  deformation  are  symaetrical  with  respect  to  the  center  0 
of  the  surface  of  contact.  Taking  Fig.  A.k  to  repreaent  the  surface  of  c«ataet, 
and  N  as  a  point  on  the  surface  of  contact  of  the  lower  ball,  the  displacement 
of  this  point  may  be  found  in  the  following  Banner.  For  a  point  load  on  the  surface 
of  an  infinite  meditsi  such  as  shown  in  Fig.  A.5  ,  the  vertical  deflection  at  a 
distance  r  from  the  load  P  Is  given  as 


w  ■ 


WBr 


(A.7) 


where  E  is  Young's  a»dulus  and  v  is  Poisson's  ratio  for  the  mediini. 
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Nov  let  UB  consider  Fig.  A. U  vhere  ve  have  a  distributed  load  q 
over  the  circular  area  of  radius  a.  Ve  are  concerned  with  the  displacement 
of  point  M  In  relation  to  the  distributed  load  q.  Ihxis  the  vertical  deflection 
V  at  any  point  M  vlthln  the  loaded  boundary  Is 


V  ■ 


(A.  8) 


vhere 


da  B  ds  .  dt  .  a 


(A.9) 


Thus 

V  -  JJ  qdtds  (A.  10) 

Hence  If  two  spherical  balls  of  equal  radii  R  are  pressed  together, 
at  SOM  point  r  the  local  displacement  v  becomes 

2w  -  a  -  r^R  -  ^  JJ qdyds  (A.ll) 

The  distribution  of  q  must  therefore  be  such  that  Eq.  (A.ll)  Is  satisfied.  It 
vill  Dov  be  shovn  that  this  requirement  Is  satisfied  by  using  a  pressux’e  distri¬ 
bution  of  q  over  the  contact  surface  represented  by  the  ordinates  of  a  hemisphere 
of  radius  a  eonstz*ucted  on  the  surface  of  contact. 

If  q^  is  the  preasure  at  the  center  0  of  the  surface  of  contact,  then 

'lo  *  a  (A.  12) 

vhere  k  Is  a  constant  factor  indicating  the  scale  of  our  representation  of  the 
pressure  distribution.  Along  a  chord  mn  the  pressure  q  varies,  as  Indicated 
In  Fig.  A.  4  by  the  dotted  semicircle. 


Perforalng  the  integration  along  the  chord  we  find 
Jqdi  .  A 


(A. 13) 


where  A  is  the  area  of  the  semicircle  indicated  by  the  dotted  line  and  is  equal 
to 


A  ■  ^  (a^  -  r^  sla^») 


Substituting  Bq.  (A.  14)  into  Eq.  (A.11)  we  find 

*  -  “  - 1'-  ^ 


d^ 


(A. 14) 


or 


a 


‘Is.  ils. 

a  E 


r^  sin^)  d* 


Integrating  Eq.  (A. 15)  yields 


a 


R  a  E 


2 

r 


) 


(A. 15) 


(a.  16) 


Ihis  equation  will  be  satisfied  for  any  ralue  of  r,  and  benc;>  ti..e  assumed  pres¬ 
sure  distribution  is  the  correct  one  if  the  following  relatl-tni;  exist  for  the 
displacement  a  and  the  radius  a  of  the  surface  of  contact  be'!^ve<!n  two  equal  radii 
spheres : 


a 


q^hjr  (1  -  v2) 
E 


q^sjr  (1  -  V) 

5* 


(A. 17) 


and 

q^FOr  (1  -  v2)  q^RT  (l  -  V) 

•  - - S -  - - STii - 


(A.  16) 


where  (li  is  the  shear  modulus  of  the  spheres 
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If  the  voluDe  of  the  pressure  dla^am  between  two  spheres  is  de¬ 


fined  as  the  nomal  force  H  between  the  spheres,  then 


R 


2 

3 


froa  which  we  obtain 


i  JL 


Coabinlng  Bq.  (A.19)  with  Eq.  (A.I7)  and  (A.I6)  yields 


The  radius  of  contact  a  may  be  ellnlnated  fron  Eq,  (A. 20)  and  (A. 21) 


to  glwe 


a 


2 


3 Cl  -  v)H 

572 


L  8MR 


-1/5 


Equation  (A. 22)  shows  that  the  relative  approach  of  the  center  of  two 
spheres  Is  a  function  of  the  two  thirds  power  of  the  contact  force. 

The  nomal  compliance  C  Is 

-1/5 


-  da  U 
‘^•di- J 


572 

Smr 


--1} 

1/2 

84R 


which  may  be  simplified  to 


C  - 


Ll--v1 

2Ma 


(A.I9) 

(A- 20) 
(A. 21) 

(A.  22) 


(A.  23) 
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where  a  i'-:  given  by 

_  [  Ml  -  v)NR 

a  -  I  gru 

L 

There  are  two  key  assumptions  in  the  Hertz  theory  which  should  be 
discussed  to  Justify  Its  use  in  a  theory  of  granular  media. 

The  first  of  these  is  that  described  by  Eqs.  (A.l)  where  it  was 
assumed  that  the  deviation  from  a  tangent  plane  is  quadratic,  i.e. 

2  . 

"  “  Ir  (A,25) 

The  exact  expression  for  this  is 


1/3 


(A. 24) 


or 


r 

2R 


1  2  1  /Z\£ 

2  R 


E 


(-)' 

'r' 


(a.  26) 


Hence  if  -  is  small  with  respect  to  1,  the  assxaaption  that  the  hi^er  order 
terms  can  be  neglected  is  quite  in  order. 

The  second  assumption  is  of  a  lower  order  than  the  first  and  re¬ 
quires  that  the  radius  of  contact  be  small  compared  to  the  sphere  radius  R. 

This  assvmiptlon  is  associated  with  the  use  of  the  expression  for  the  deflection 
resulting  from  a  point  load  on  an  infinite  medium,  I.e.  Equation  (A. 7).  Eesentlally 
>rtiat  is  being  assvined  is  that  the  curvature  of  the  spheres  beyond  the  contact 
area  does  not  affect  the  stress  distributions  and  deflections  at  the  contsicts. 
Schematically  this  assvnptlon  implies  that  the  surfaces  In  contact  are  two 
infinite  media  with  small  bumps  In  them  rather  than  two  sphere  (Fig.  A. ^).  Clearly 
if  the  radius  of  contact  is  small,  the  curvature  has  little  effect  on  the  behavior. 
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This  Is  partlculeurly  tnte  vhen  the  tw>  contact  surfaces  eire  spherical  and  the 
stresses  avay  froa  the  contact  area  damp  out  very  fast.  Ihe  tvo  problems 
Illustrated  by  (Fig.  A. 3)  could  conceivably  be  solved  and  a  quantitative 
comparison  be  made;  however,  the  effox^  would  be  a  major  one  And  the  value  to 
the  study  at  hand  Is  questionable  as  Is  seen  In  the  following  discussion. 

Ihe  problem  of  concern  has  to  do  with  the  use  of  the  Hertz  theory 
to  describe  the  behavior  of  a  granular  medium,  namely  sand,  subjected  to  an 
average  applied  stress,  0^,  over  the  surface,  therefore.  It  Is  of  Interest 
to  determine  vbat  limits  are  to  be  placed  on  this  stress  due  to  other  reasons 
such  as  crushing.  In  fact  It  Is  known  that  cmshlng  begins  In  sand  compressed 
one  dimensionally  at  average  stresses  on  the  order  of  1000  psl.  Hy  the  time 
the  stress  reaches  3,000  psl  the  crushing  effect  Is  a  major  part  of  the  deformation 
behavior. 

In  Section  3'1'1  it  \i3.s  determined  that  the  normal  contact  force 
7<2  '^■^‘’•■2  related  to  the  radius  of  contact  by  Eq.  (3.I) 

i 

^  (A.27) 

The  normal  force,  Ii2,  was  related  to  the  applied  force,  P^,  by 
Eq,  (3.48a) 


3  (1  -  v)  RN^ 
8^1. 


No  (1  + 


i;  -  g 


Pz 


(A. 28) 


Furthermore  it  vas  shovra 
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that  the  average  stress  a  Is  related  to  P  by 

z  z 


0  » 
z 


SR-i 


(A. 29) 


Equations  (A. 27),  (A. 28)  and  (A. 29)  may  be  combined  to  give 


(-) 

'R'' 


2^^(l  -  V  )  ^ 

8(1  +  f) 


(A.  30) 


Ibe  relationship  between  a/R  and  the  stress  Is  Independent  of  R. 
Representative  values  of  f,  v  and  ^  for  qxjartz  which  are  appli¬ 
cable  to  sand  cure 

f  =  0.15 
V  =  0.20 
U  =  6  X  10^  psl 

Using  these  values  and  Eq.  (A. 30)  we  obtain  the  following  re¬ 
lations  between  a/R  and  a  : 

2 


0.369  ^ 


(A.31) 


a 

R 

a 

2 

pel 

O.OJl 

500 

0.039 

1,000 

0.057 

3,000 

0.085 

10,000 

0.10 

16,500 

Thus  we  can  see  that  for  the  stresses  of  Interest  (prior  to 
crushing)  where  the  theory  is  to  be  applied,  the  radius  of  contact  com- 
pcured  to  the  radius  of  the  sphere  will  be  much  less  than  0.10.  Rot 
only  is  this  fairly  gman  tut  it  should  be  remembered  that  there  are 
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Other  assumptions  involved  in  the  application  of  this  theory  to  sand 
vhich  may  be  far  worse.  In  fact,  the  assumption  that  sand  particles, 
even  uniformly  smooth  ones,  are  spheres  is  perhaps  more  questionable 
than  this  one.  Microscopic  photographs  of  sand  particles  show  many 
types  of  asperities  irtilch  naturally  cannot  be  easily  considered  in  this 
theory.  If  the  effort  warrants  it,  it  may  be  possible  in  the  future  to 
extend  the  theory  to  include  rouj^ess  effects  by  using  the  recent  work 
of  Goodman  [A.3].  Such  an  extension  is  not  warranted  at  the  present  in 
view  of  the  nonuniformity  of  sand. 

Another  effect  ignored  is  that  of  the  tangential  forces  at  the 
contacts  on  the  geometry  Just  outside  the  contact  area.  With  large  con¬ 
tact  forces  there  may  exist  little  outward  bumps  outside  the  contact  sur¬ 
faces  which  would  have  some  effect  on  sliding. 

From  the  above-mentioned,  it  can  be  concluded  that  within  the 
confines  of  the  expected  use  of  this  theory,  the  assumption  the,t  the 
Hertz  theory  is  applicable  to  the  behavior  of  granular  materials  is  cer¬ 
tainly  warranted. 
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Fig.  A.  1  TV/O  UNSTRESSED  SI'HERES  IN  CONTACT 


Fig.  A.  2  SEGMENT  OF  A  SPHERE  WITH  A  TANGENT  AT  0 


Fig.  A.  3  POINT  LOAD  ON  AN  INFINITE  SURFACE 
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APPENDIX  B 


NUMERICAL  SOLUTION  TO  THE  EQUATIONS  GOVERNING 
THE  BEHAVIOR  OP  A  GRANULAR  MEDIUM  HAVING 
RESTRICTED  LATERAL  STRAINS 


Ihe  general  equations  governing  the  Behavior  of  a  face  centered 
array  of  spheres  subjected  to  a  vertical  force  and  symmetrical  lateral 
forces  P„  have  been  developed  in  Section  3-l"l"  "nie  resulting  equations 

n 

including  the  lateral  strain  condition  were  given  previously  as  Eqs.  (5.37) 
and  take  the  following  form: 


dNg  +  dTg 


/2 


sTs 


dN^  +  dNg  -  dT^  -  dPg  =  0 


^  dN-  -  dN  +  ^2  dT_ 
^2  ^2 


dN, 


He 


{B.l) 


Here  the  N's  and  T^  are  the  normal  and  tangential  contact  forces 
and  the  C's  and  S^  are  the  corresponding  compliances.  (Pigs.  B.l  and  B.2). 
The  radii  of  the  si^eres  is  R  and  the  lateral  strain  is  c^.  Let  us  intro¬ 
duce  the  dimensionless  variables 


^  R^,  ^ 


3(l-v)  ”2 


?(l-0 

8 


Ph 


8 


(B.2) 


where  ^  and  v  are  the  shefu:  modulus  and  Poisson's  ratio  for  the  spheres. 


1S6 


Substltutljag  Equations  (B.2)  into  Equation  (B.l)  yields 


dttg  +  dtg 


-  T*®* 


/2 


+  dng  -  dtg - j  dp. 


H 


2-V 


^  dn^  -  dng  +  gQ-.vT'  “^2  “  ° 
“l 


3  3/5 
2“l 


dn. 


idiere  K  depends  on  t.  and  as  foUovs : 

£  Otg 

a.  ?0r  tg  Increasing 


(B.3a) 

(B.3b) 

(B.5c) 

(B.Jd) 


*“2  1 


dn-  ,  dn-  dn_  t-  -l/3 


b.  Eor  t  decreasing 

2 


*^2  1 


dn^ 


dn. 


dn. 


t*  -  t,  -1/3 


2ftj- 


S) 


(B.5) 


Idiere  t«  is  the  hipest  ralue  obtained  bgr  tg  and  is  the  Talus  ftroa  idiieh  it  Is 
decreasing. 


Bie  strain  in  the  Tertieal  or  s  direction  is  glTea  bgr 
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dc. 


do^ 


2-v 

5(1^ 


Kdt. 


(B,6) 


•ad  the  stressee  la  the  rertleel  and  horisontel  dlx«etloiia  ere  glrea,  reepeetlrely, 

By 


®H"3Tr7T 


(B.7) 


Ihe  problee  ie  nov  to  detexalne  a  eolvtioa  to  the  Bqmtlons(B.3), 

Iheee  eqaatloae  are  dirterential  e^pwtioaa  la  texae  of  the  inflniteelaal 
laerieate  of  the  forees.  One  aethod  of  eolvtioa  la  to  latesrate  theae  e^vatloaa 
exactly  ahere  each  aa  late|^tte  eaa  be  done.  Sveh  aa  iatagration  la  carried 
out  la  Section  for  the  caae  ahere  the  lateral  atrala  la  saro.  Qafortnnately 
It  does  not  appear  that  aa  exact  aolutlon  caa  be  obtained  for  other  eaaea, 
particularly  for  the  cycle.  Bence,  la  order  to  atudy  the  btiiarlor  of  the 

nedlm  acne  other  nethod  of  aolvtloa  naat  be  carried  out. 

Iquatlona  (B.3)  are  Ideally  auited  for  anierleal  aolvtloa  on  a 
conpnter.  If  ac  coaalder  the  equatlOB  not  aa  dlffareatlal  laeranenta  d  but  aa 
finite  laemaaata  A  the  eguatlona  beccne 

hn^  +  Zing  -  0*2  -  -  0  (B.8) 


-dOg  +  lgiTyBOtg-O 
“l 

5  V3 

^  -fn^  0«H 

ahere  K  la  glTea  by  Iquatlona  (B.h)  or  (B.3)  vlth  dlffereneaa  raplaclng  the  dlf- 


ferentlala 
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If  the  initial  conditions  are  knom^  Bqtiations  (B.8)  can  he  used 
to  trace  the  behavior  of  the  spheres  through  a  CTcle  of  loading.  Hie  Tertieal 
load,  p^,  the  lateral  strain  or  the  lateral  load  Pg  aagr  be  increaented  end 
a  solution  to  Equations  (B.8)  gires  the  inereaents  vhieh  occur  to  the  other 
forces  as  a  restdt  of  these  changes.  Likevise  froa  these  results  the  inereaent 
to  the  rertical  strain  aagr  be  detexained  Area  Equation  (B.6). 

It  should  be  noted  that  Equations  (B.8)  are  non>linear  difference 
equations  because  K  is  a  fuaetion  of  n^  and  t^.  Bence,  the  equations  Bust  be 
solTed  by  iteration  for  each  load  inereaent  cycle. 

It  was  also  noted  in  Section  ^.2  that  these  equations  hare  a  singularity 
at  sero.  therefore  it  is  not  possible  to  obtain  any  rational  b^utnrlor  by 
starting  with  the  noraal  forces  at  the  contacts  equal  to  sero.  TO  offset  this 
problea  «e  introduce  an  Initial  hydrostatic  stress  p^  vhieh  in  turn  creates 
an  initial  noraal  force  at  the  contacts  of  n^.  As  vas  shown  in  Section  ^.1.2  this 
gives  stress -strain  and  contact  force  relationships  as  follows: 


O  0 


(B.9) 


(B.10) 


where 


®o"  jTT^o 
3(l-v)  p^ 


(B.11) 


3(l-v)  S 


(B.12) 


The  addition  of  this  Initial  stress  condition  does  not  alter  the 
governing  Eqxaations  (B.8)  pxwided  that  the  total  nomal  contact  forces 
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used  in  the  eonpllenee  texas  and  strain  equations  ineludle  the  hydrostatic 
effect.  Hence  all  variables  except  the  tangential  forces  include  the 
hydrostatic  contribution.  Ihe  behavior  in  addition  to  the  initial  hydrostatic 
effect  or  the  total  effect  can  then  be  obtained. 

In  the  eoBputation  procedure  we  remsber  that  the  theory  is  restricted 
to  that  bdunrlor  where  sliding  does  not  occur  at  the  contacts.  Ibis  is  re¬ 
flected  in  the  condition  that  K  aust  be  greater  than  sero.  VAien  K  becoaes 
negative  (or  zero)  sliding  occurs  and  the  solution  is  no  longer  valid. 

Furthezaore,  to  speed  convergence  of  the  iteration  process  and  to 
take  care  of  any  divergence  during  the  ccaputatlon  for  one  inereaent,  a  con¬ 
vergence  subroutine  was  included  in  the  prograa. 

A  step  by  step  solution  is  as  follows: 

(a)  Increaant  loads  and/or  strains 

(b)  Calculate  all  constants 

(c)  Calculate  dn^  froa  Equation  (B.Jd) 

(d)  Assrae  values  tor  K  and  ^h/n^  on  the  basis  of  the  previous  cycle 

(e)  Coapute  dt^  froa  a  co^ination  of  Equation  (B.^)  and  (B.^c) 

it)  Deteraine  dn^  froa  Equation  (B.}a) 

dn- 

(g)  Check  T*  for  greater  than  or  less  than  l/f 

dtg 

(h)  Coamte  K  froa  Equations  (B.h)  or  (B.^)  depending  on  the  results 

of  (g).  If  K  is  negative  then  stop  since  the  equations  are  invalid. 

(l)  Coiqpute  — 

(j)  If  assuDWd  in  (d),  add  increaents  1 1  total 

variables  and  go  on  to  next  cycle 

(k)  If  “l/ug  do  not  agree  with  assused,  repeat  steps  (f)  thro- 

(l)  using  results  of  this  cycle. 
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A  flov  dlagraa  of  th«  cooputation  sequence  Is  given  as  Fig.  B.3« 

In  order  to  study  tbe  accuracy  of  the  eaqmter  solution  and  to  compare 
Its  results  vltb  a  closed  foxm  solution,  a  solution  vas  obtained  to  the  one 
dimensional  problem  (c^  ■  O).  ihls  corresponds  to  that  developed  In  Section 
3.2.1  for  idileb  an  exact  solution  has  been  obtained,  namely.  Equation  (3«3l). 

Figure  B.4  shove  a  comparison  of  the  results  obtained  from  the  computer 
program  vlth  those  resiLLtlng  from  Equation  (3.3l).  accuracy  vas  excellent 

and  DO  difference  Is  discernible  on  a  graphical  plot.  The  closeness  of  the  tw> 
results  la  Illustrated  by  the  sample  values  noted  on  the  figure.  Further 
vork  Is  anticipated  on  this  phase  of  the  program  and  vlU  be  carried  out  as  the 
need  arises.  It  la  expected  that  It  vlU  be  useful  in  studying  the  unloading 
phases  and  those  eases  vhere  some  lateral  strain  Is  alloved. 


I 
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Fig.  B.  1  UNIT  CUBE  OF  A  FACE  CENTERED  CUBIC  ARRAY 

OF  EQUAL  SPHERES  SUBJECTED  TO  INCREMENTAL 
FORCES  IN  ONE  DIMENSIONAL  COMPRESSION 


Fig.  B.  3  FLOW  DIAGRAM  OF  THE  COMPUTER  PROGRAM 
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